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High order terms in the electromagnetic multi-pole development expose a stabilizing 
mechanism for the atomic orbitals in the presence of a random background of elec- 
tromagnetic fluctuations. Boyer and Puthoff set forward the idea that for the Bohr orbits in 
the hydrogen atom, radiation losses could be compensated by absorption from the QED-predicted 
Zero-Point Field (ZPF) background. This balance is, on average over the orbit, a necessary condi- 
tion for stationarity of the orbits, imposing a relation on the pair Ro (orbital radius), uio (orbital 
(N- angular velocity) ; such relation is simply what we have for long known as angular momentum (AM) 

quantization [I — 1). Nothing has been said yet, however, on how could this balance be attained 
on a quasi instantaneous basis, in other words, how could the orbit accommodate the instantaneous 
excess or defect of energy so as to keep constant the (at least average) values of its parameters (Ro, 
ujo)- Using classical electromagnetism, we explore some high order interactions between realistic 
particles, exposing a mechanism (a feedback loop between variables) that makes that stability pos- 
sible. Puthoff 's work led necessarily to the quantization of AM: "if stable orbits exist... 
then their AM must be quantized"; now we are able to do a much stronger statement: 
"the equations of the system, in the presence of ZPF background, lead necessarily to a 
O | discrete set of stable orbits". 

The job is done in two steps. First, we work with an ideal electromechanical model composed of 
two (charged) particles: one is a point-like magnetic source and the other one consists of a spatially 
"extended" charge distribution, orbiting around the point like one. In the second place, a relocation 
of the inertial reference system and other minor changes are done to reinterpret the former picture 
in terms of a realistic hydrogen atom, the two situations turning out to be formally identical. As a 
difference with traditional approaches, here it is the inner structure of the nucleus ( and 
not of the electron) that plays the crucial role in our model. The following step regards the 
existence (for each I) of an infinite and discrete spectrum: the presence of a secondary feedback 

■ loop in the equations is crucial for the feasibility of excited states. Up to the order of our 
equations and under some further approximations, our model admits a continuum of stationary 
trajectories described by four parameters: {Ro,oJo, Ri,i^i}', the first pair (Ro,Uo) corresponding to 

^s^j . the orbital movement, the latter (Ri,u>i) to what we have called a "secondary oscillation" (may 

' there from here a connection to Cavalleri and others' model of spin as a residual helical oscillation 

of a point-like electron [2S|?). Stationarity determines, for a given Ro, both ojq and Wi (imaginary 

■ poles in the linearized frequency description of the system), while Ri remains free, allowing, once 
radiation/absorption (rad/abs) processes are considered, the power balance that is necessary for 
stability, while discreteness of the spectrum is now retained via an additional condition, a phase 
relation Wi = n/uio (n integer) that we regard also necessary for stability. Such phase relation also 
allows us to arrive to an E-spectrum that could be made to resemble, under some approximations, 

1 the well known 1/n 2 one. Of course stationarity is so far only a necessary condition for stability: 

more as an intermediate step than a rigorous proof, we introduce a stochastic description of 
the rad/abs processes, in a way that makes (at least) possible the negativity of the real 
part of the eigenvalues (poles in the frequency domain) of the linearized description of 
the system. 

Obviously, circular trajectories can only give rise to p-orbitals (or at least I > 1 ones): an 
additional (infinite and discrete) set of trajectories with no net AM can also be found, 
that would account for the s -spectrum. Both in the s or p cases, an atomic transition implies 
an emission/ absorption of a wave-packet with energy h(u>i —uij), being the (main) 

frequencies of the oscillation (primary osc.) of the orbitals involved. Besides, for the 
feasibility of our model, we are led to think in fully relativistic orbitals, "nodes" absent at least 
for / > [13]; anyway, the unrealistic nature of pure states would also solve this problem, even for 
I — (Fritsche). Finally, the action of the stochastic background on the oversimplified, entirely 
deterministic orbits that we provide here would produce a probability distribution extended to the 
whole space, an estimation which is left for elsewhere. We also barely touch other issues like the 
extension of the results to 3D, as well as to / > 1 orbits; providing a reasonable framework for future 
developments is, in any case, the ultimate goal of this paper. 
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Previous note: QM and local theories 

In the past we have come across real difficulties to make 
anyone, not only read our work here, but even just discuss 
about the main idea behind it as a purely electromechan- 
ical problem. Such idea, whatever inaccuracies or details 
we may have overlooked in the following pages, works 
in the abstract level: including additional, "inner" de- 
grees of freedom in the description of a classical system 
is an absolutely necessary step to make stable behav- 
ior possible, as whoever is acquainted with the basics of 
Dynamical Systems Theory will not find difficult to rec- 
ognize. We were warned by Dr. Risco-Delgado about 
this long ago, but took us nevertheless quite longer to 
realize the gravity of the problem: we had that "QM is 
non-local" barrier that relegated to utter marginality any 
effort aiming at, or simply suggesting, the possibility of 
reconciliation for the classical and quantum worlds. 

Well, recently we have make quite an effort, our own, 
to examine the question in detail. Our conclusion: QM 
may be non-local; nature is not, at least based on the 
experimental evidence of the last decades, and of course 
until the opposite is proven (and at this point it really 
does look highly unlikely [39j). According to us, then, 
there is no such unsavable boundary between the 
two worlds, time to start giving some credit to other 
alternatives: they may better what we have (and which 
works) or they may not, but in any case they will help us 
understand with more profundity... is it not what science 
is about? 

But that is not only our own conclusion: aside from our 
own recent advances in what we call "the Wigner-PDC 
picture of photon entanglement" [37], IH[ , we suggest the 
reader take a look at the last developments by De La 
Peha and collaborators, a paper named "Origin and 
meaning of non-locality" (2011). According to us [36| . 
and as far as we can see, also in principle according to 
them, Quantum Mechanics is not a theory on its own 
right: QM is a shortcut for a much more complex, 
entirely classical world of particles and fields. 
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I. PROLOGUE 

The aim of this paper is to show that an elementary 
mechanism exists, in the context of SED (classical elec- 
tromagnetism plus a zero-point background of radiation), 
that makes not only possible but necessary the discrete 
and stable character of the hydrogen orbital spectrum, 
and by extension the atomic one in general. Such a 
stabilizing mechanism arises in a higher order descrip- 
tion of the system, where an elementary model of the 
nucleus with some inner structure is included. When 
complexity is present in the structure of a systems, it 
is not a rare thing that its dynamics presents some sta- 
bility and/or attraction phenomena, and very often far 
beyond the obvious. Therefore, beside customary classi- 
cal electromagnetism and the interaction with a stochas- 
tic background, the inclusion of an elementary model of 
the nucleus stands for that new ingredient in the pic- 
ture: "complexity" . To advance some evidence on the 
feasibility of this approach, we have the proof that our 
picture can accommodate naturally (something close to) 
the well known ^-dependence of the hydrogen spectrum, 
though this is done on the basis of a somehow arbitrary 
frequency matching condition that we have not been able 
to justify convincingly, at least so far. The space for a 
feasible explanation is there, however. 

Besides, the material this paper attempts to cover is 
quite extensive. We think therefore that the most conve- 
nient way to start is to provide a brief, but more or less 
complete overview of the whole framework we are trying 
to set here. This section can be considered, in a way, as 
an "extended abstract" , making the exception that some 
of the ideas that we mention now are not yet sufficiently 
formally developed or may even not appear again in the 
rest of the document (although we think they are suffi- 
ciently relevant to include at least a mention) . This said, 
we will be quite satisfied if our main results are under- 
stood at least to the point that they can be a subject of 
debate. 



A. Settling some foundations 

The following ideas are, somehow, both a point of de- 
parture and arrival: 

(I) To be able to relate the concepts of quantum and 
classical angular momentum (AM), we will assume that 
each projection of quantum (orbital) AM corresponds to 
an average on (the projection of) the standard classical 
AM along a closed, periodic trajectory. This is justified 
in Sec. IXVI At it seems almost obvious to us given the 
corresponding addition rules in QM. 

(II) On the other hand, for us point (I) also applies to 
the inner AM or spin of a particle. Indeed, we will often 
talk of a "classical spin" when referring to a rotational 
movement of a particle (or a distribution of particles) 
around an axis of symmetry. Whether if quantum spin 
appears in a representation where an isomorphism can 



be established (or not) to the group SO(3) of rotations 
in ordinary space is a question that we have analyzed 
somewhere else [4(3]; however, here we are only interested 
in the fact that, whenever a particle has spin, at least 
some sort of classical AM can be associated to it (and in 
particular in regard to its magnitude). Besides, we are 
also aware there are some relativistic difficulties with the 
picture of a classically spinning sphere of charge; we do 
not consider them definitive though [33| , nor is "rigidity" 
of such sphere nothing more than merely a simplifying 
assumption. 

(III) In relation to the former two points, we will be 
thinking in QM as a "semi-static" theory that masks a 
richer dynamics underneath, involving perhaps a higher 
number of (hidden) classical degrees of freedom (that 
"complexity"): for instance, the quark dynamics (as 
point-like entities carrying charge and spin, but configur- 
ing a distribution, with an associated dynamics, in space) 
inside a nucleon. 

(IV) In coherence with (III), while average values of 
the projections of those completely classical AM will stay 
attached to a discrete spectrum as in QM, their instanta- 
neous values (and in a micro-dynamics that is transpar- 
ent to atomic transitions) can, at least, oscillate around 
the quantum mechanical one. Interaction with a random 
background could cause much of those fluctuations, and, 
incidentally, spontaneous state transitions that are noth- 
ing new in orthodox QM. For instance (and we may be 
going unnecessarily too far, too soon), the orbital AM of 
the quarks inside a proton (therefore adding to the parti- 
cle's inner AM or spin) can oscillate around its quantum 
value, in response to perturbations coming from the inter- 
action with the background. Because those oscillations 
are more or less minimal, we propose the term "residual 
spin" (RS). 

(V) The appearance of the "photon" (E — hut) as a 
constraint on the spectrum of these possible "discrete" 
exchanges of energy-momentum between metastable 
states of the systems... seems completely natural in this 
framework. 

With all this in mind, we insist that will be dealing 
here, exclusively, with classical Maxwellian electromag- 
netism. We have charges (or distributions of them) and 
fields, amongst them a random background: fluctuations 
of the value of the fields within the vacuum. The ref- 
erence to the quark model is necessary (is solves some 
difficulties of the model with special relativity), but nev- 
ertheless nothing more than tangential. 



B. Warning on names and conventions 

Besides, we will need to adopt, unavoidably, some new 
names and conventions, that may undergo changes in fur- 
ther versions of this paper. In particular, we are far from 
satisfied with the term Inverse Magnitude Spin Or- 
bit coupling (IMSO). Amongst former "candidates" 
was the term "classical magnitude-driven (proton) spin 
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- (electron) orbit coupling". Not being completely con- 
vinced and for the lack of a better choice, we have de- 
cided, for the moment, to stick to something shorter. It 
happens so with the term "(proton) residual spin", that 
we use to refer to a residual inner AM coming from the 
orbital movement of quarks inside the proton. In other 
words, by this term we refer to a classical rotational 
movement of the proton "as a whole" , around its axis 
of symmetry. The need to consider this kind of motion 
will become clear (we hope) later. Indeed, it is the "oscil- 
lations" in this movement that will make our stabilizing 
mechanism work. 

Besides, perhaps misleadingly we have also indulged in 
trying to present our new terms in the multipole devel- 
opment of the hydrogen as related with the well known 
LS term in atomic physics (and its classical counterpart). 
Maybe we were influenced by the prominent role played 
by that term in relation to stability in Nuclear Physics, 
and we also had the intention of providing the reader 
with a more familiar context. Although that relation ex- 
ists, they in any case terms of a very different nature, 
a recollection of differences made, for instance, in Sec. 
IXVIII Bl Some naming decisions have been made in that 
direction, too, which now seem not so appropriate. 



C. An overview of results and some general context 

We provide a very brief summary of the paper. Basi- 
cally, what we do is: we begin this effort by analyzing an 
elementary electromechanical model. Later, we make the 
necessary changes to apply our results to an elementary 
model of a hydrogen atom. 

Under the action of an ideal magnetic dipole, a moving 
charge with some spatial extension (non- vanishing second 
order moment) experiences a force and a torque, medi- 
ating between its "orbital" and "self rotation" degrees of 
freedom, this last what we will call a (classical) "residual 
spin" (RS). That pair (force, torque) introduces a bidi- 
rectional coupling between the (instantaneous) values of 
two classical angular momenta corresponding to those 
two degrees of freedom, a coupling that acts on their 
magnitude, as a difference, for instance, with the spin- 
orbit (LS) coupling in atomic physics, which acts only on 
their relative orientations: this is no surprise but only a 
consequence of the fact that our mechanism corresponds 
to terms of a higher order in the multipole development. 
For instance, too, the LS coupling only involves a punc- 
tual value of the field, therefore being unable to feel its 
gradient, which is another fundamental difference with 
the interaction we expose here. As we said before, we 
have named it, provisionally and for the lack of a bet- 
ter choice, Inverse Magnitude Spin Orbit coupling 
(IMSO). 

An interpretation of our model in terms of the hydro- 
gen atom is more elaborate, but we conclude that similar 
results could be applicable to the orbitals of the electron 
around the nucleus, when a spatially "extended" charge 



distribution is associated to this last, and always in the 
presence of a background. On the other hand, later, to 
be able to interpret our result in a more realistic situation 
(an atomic model), the role played here by that classical 
spin will be taken by what we have called residual spin 
(RS) of the nucleus: in the simplest case of a proton, 
a residual (inner) AM coming from orbital movement of 
quarks inside its structure. 

In the presence of a certain background compensating 
on average the radiative losses, we are already able to 
prove that the identified "feedback loop" (FL) between 
orbit and the RS-N provides a necessary and sufficient 
condition for the existence of stable "orbitals" in the clas- 
sical configuration space of the system. These orbitals 
are produced by oscillations around a set of privileged 
trajectories or "attractors" , that (we have proven) can 
either bear or not net (average) AM. For instance, this 
last would be the case of the s-orbitals, with vanishing 
average AM. Both the stability of atomic orbitals and the 
discrete character of the spectrum can be, then, at least 
potentially explained. 

Curiously enough, our preliminary calculations show 
that the dependence of the associated energy correction 
in the (classical) orbit radius agrees with that L ■ S term 
in atomic physics. Nevertheless, those energy contribu- 
tions are not net ones: they are only translation of energy 
from one degree of freedom to another, and moreover, 
their average (seems to) vanish on a whole cyclic tra- 
jectory (an orbit). This is related to the fact that all 
new terms arise from the Lorentz law, whose associated 
(non-conservative) work on a particle is always zero. As 
a result of this, no new observable contribution to the 
spectra appears. 



D. Relation with the Bohr model 

We depart from the assumption that quantum angular 
momentum (AM) is related to an average classical AM, 
along a closed trajectory. This is justified in Sec. IXVI Al 
Aside from other differences, our approach veers from the 
Bohr approach in two important points: 

(i) The Bohr approach works with the purely Coulomb 
potential. 

(ii) The Bohr approach only seeks for circular or ellip- 
tical orbits. Clearly, this cannot account for the actual 
ground state of the real hydrogen (an s-orbital in QM), 
because Bohr's trajectories show net AM. In our frame- 
work, we do find a set of stable trajectories with van- 
ishing (average) AM. Those trajectories, including both 
the ground state, introduced in Sec. [Xj and its excited 
states, see Sec. IXI1 are not possible in a problem with 
a purely electrostatic potential. We have to resort here 
to terms related to the Lorentz force acting on an "ex- 
tended" model of a particle. Details on all this will be 
given in the following pages. 
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E. Other additional questions 

To complete this momentarily vague picture, at least 
two other questions should also be mentioned: 

(I) Relativistic considerations: the need to consider a 
(classical) "residual" spin (RS). We are aware of the dif- 
ficulties that arise from assuming, for the proton, the 
model of a uniform sphere of charge, rotating with an- 
gular velocity proportional to its phenomcnological mag- 
netic moment. These difficulties, of relativistic origin, 
can be safely ignored if a more realistic model of quarks 
(carrying both charge and magnetic moment) is assumed. 

If the quark model is assumed, the contribution of the 
angular velocity of the proton as a whole, around its axis 
of symmetry, to its overall inner AM is either negligible 
or at least vanishes on average. However, classically, the 
magnitude of this magnetic moment can still oscillate. 
We have renamed this oscillation dynamics as "residual 
spin" (RS), and it can play the role that, in our first ide- 
alized situation, conventional (classical) spinning move- 
ment played itself. This RS would eventually correspond 
to oscillations of the orbital AM of the quarks inside the 
proton (oscillations with no quantum counterpart, but 
with complete sense in our classical reformulation of the 
system) . 

Besides, we make several references to the "rigidity" 
of this distribution, but this is simply a device to make 
our argument clearer. Indeed, the mechanism we expose 
is dominant up to a certain order, whenever the two ob- 
jects remain sufficiently apart. Higher orders may deform 
the shape of the distribution, or have other effects that, 
whenever the range of distance is the appropriate one, do 
not have to bother us. 

(II) The 3D problem. We will provide some preliminary 
ideas about the extension of the argument to three di- 
mensions (3D). In particular, that extension would need 
to include precisely the (also well known) classical coun- 
terpart of that quantum LS interaction, acting as a mod- 
ulator of the strength of our "stabilizing mechanism" 
(IMSO). On a leading order, their combined effect would 
allow for a "modulated precession" of the axis of the clas- 
sical orbit, necessary to account for / > orbitals (there- 
fore associated to a certain privileged direction in space). 

F. But why would this be interesting? 

Two kinds of reasons. The first, scientific hygiene. 
This is not an eccentric theory with all kinds of new enti- 
ties, dimensions and laws of interactions with an obscure 
meaning. This is a set of charges and Maxwell's laws. 
Someone has to see how far it (really) can go. Someone 
has to do it, but, yes, we admit this alone may make it 
convenient, but hardly interesting if it does at all. The 
second, the possibility of going beyond the standard set 
of phenomena that QM describes so economically well. 
Against skepticism, here are some hints that may lead to 
somewhere: 



(a) A model theory built to account for the mechanism 
we expose (for instance based on an effective potential) 
would not suffer from renormalization problems: the kind 
of mechanism that we describe here constitutes a "leading 
order" behavior. In fact, it gets smeared out by higher or- 
der terms when the point particle and the distribution of 
charges get too close. When too far apart, naturally, the 
interaction also becomes negligible. This kind of behav- 
ior is very suggestive, as imposing a lower limit on the 
distances is equivalent to introducing the conventional 
cut-off in momentum transfer which is a common rule 
(trick) in Quantum Electrodynamics (QED). Therefore, 
if we reformulated our mechanism in terms of a new "ef- 
fective" interaction, it would not need to be renormalized 
at any scale (as we said, through the inclusion of an en- 
ergy cut-off, for example). 

(b) Moreover, and though we do not plan to go any- 
where near for the moment, we are led to think: could we 
get rid of all renormalization problems in QED (and clas- 
sical electromagnetism) if we regarded it as an effective 
description of a simplest interaction law, but involving 
more degrees of freedom? For instance, could we gener- 
ate electromagnetism (with its vectorial potential) from a 
bare scalar potential acting on particles with some extra 
inner degrees of freedom (for instance, instead of being 
point-like entities, let us allow them to deform elasti- 
cally)? 

(c) Last but not least: could our "secondary oscilla- 
tion" be related to Cavalleri and others ' model {291] of the 
electronic spin as a helical oscillation of a structureless 
electron around its main trajectory? 



II. ATOMIC STABILITY AND SED 

In this section, we provide some academical back- 
ground where to set our work, as well as describe the 
content of the rest of this paper. 



A. SED: Stochastic Electrodynamics 

Our account of SED is not very accurate in regard to 
chronology and contains important omissions: an 
exhaustive review can be found in WH l. 

The electromagnetic radiation of an accelerated charge 
predicted by Maxwell's laws is often cited as the (obvi- 
ous) reason why the hydrogen atom cannot be classically 
stable, and one should then invoke the quantum mechani- 
cal Pauli principle to account for the existence of ordinary 
matter, at its common atomic phase. In spite of this fact, 
the reader will surely agree that at least the question of 
how unstable the atom is, or, how quasi-stable it may be, 
from strictly classical grounds, still retains interest. 

Indeed, much in the spirit of the, already old, theory 
of Stochastic Electrodynamics (SED), started by Braf- 
fort, Marshall and Boyer (amongst others), several at- 
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tempts have been made to explore this question in de- 
tail. SED's main addition to the classical picture of 
electromagnetism is the presence of a random, homoge- 
neous, isotropic and Lorentz invariant electromagnetic 
radiation that permeates all space even at a zero tem- 
perature (refs), the so-called zero-point field (ZPF). This 
ZPF background coincides with QED predictions for the 
vacuum. 

Marshall Q and Boyer Q , based solely on that Lorentz 
invariance, established the spectrum (with a cubic depen- 
dence on frequency) of this ZPF, and, as early as 1969, 
T. Boyer was able to derive the famous black body ra- 
diation law of Planck, without the need to assume any 
quantization Q. This very recommendable paper has 
been followed, over the years, by some other very sugges- 
tive ones, mainly by Marshall and Boyer themselves, and 
D.C. Cole. For instance, we are aware of and 
also of [rZrll. 

While, on theoretical grounds, the very origin of the 
ZPF was also investigated by Puthoff in Q , there is also 
some evidence on the experimental ones. Specifically, 
the observed Casimir force between two neutral metal 
plates in the vacuum is itself evidence for the existence 
of some kind of background. In any case, as we said, 
the existence of ZPF constitutes the main hypothesis of 
what is come to be known as Stochastic Electrodynamics 
(SED). We will also refer to [2l| for a brief account of 
SED "successes". 



B. SED (alone) is not enough 

However, up to now, these ideas have not been enough 
to completely build a convincing bridge from the classical 
to the quantum theory, as serious difficulties appear with 
a broad range of quantum phenomena. The stability of 
atomic orbitals is one of them: in this context, recent 
attempts have been made to examine in detail the inter- 
action of a simple system, such as the hydrogen atom, 
with certain kinds of radiation. Some very interesting 
attempts use a numerical approach |22h26| . Although 
an impressing piece of work, they cannot be considered 
anything further than preliminary explorations. That is 
also the stage, a preliminary one, where we are content 
to place this work. 

From our point of view, a particularly relevant effort is 
the one by Puthoff [2l[ , who explains the ground state en- 
ergy of the hydrogen atom from a dynamic equilibrium 
between the radiation emitted by the accelerated elec- 
tron and the radiation absorbed from the zero-point field 
(ZPF) fluctuations of this electromagnetic background. 
It is particularly remarkable that discreteness of the AM 
spectrum arises here as a natural and almost obvious con- 
sequence. However, in spite of that and the fact that it 
avoids, on average over time, the radiative collapse of the 
Bohr atom, it does not explain how that equilibrium may 
exist on a quasi-instantaneous basis. It does not explain, 
either, why there should exist a discrete set of energies 



for the orbitals with a given AM. To this matter we aim 
most of our efforts here. 

Before we continue, we feel the need to insist once 
more on the fact that we will be using here just classical 
Maxwellian (hence fully relativistic) electromagnetism, 
and, really, nothing more. Indeed, the picture of SED 
described above comes here simply as an inspiration to 
us. The kind of mechanism we propose retains its in- 
terest (and may work, too) regardless of what particular 
spectrum we choose for our radiation field, the only el- 
ement needed is the capacity of a system to exchange 
power with a certain background, both for dissipation 
and absorption. 

Yes, an accelerated electron radiates, but also absorbs 
radiation (as already pointed out by Puthoff). 

According to SED, this could be the ZPF, but it is not 
at all a necessary assumption in the rest of this work: 
any kind of radiation, for instance electromagnetic noise 
coming from the rest of the charges in the universe, would 
work. Therefore, in principle the loss of energy could be 
balanced by the absorption of power by contact with a 
bath of radiation. As we said, the idea is not new, and 
others have tried to study and simulate the behavior of 
a system in limited ranges of this kind of circumstances. 
We find, indeed, a very graphical description of results 
in [26| : "a detailed simulation of the effects of classical 
electromagnetic radiation acting on a classical electron 
in a classical hydrogen potential, results in a stochastic- 
like motion that yields a probability distribution over time 
that appears extremely close to the ground state probabil- 
ity distribution for hydrogen" . 

But a bath of radiation, whether or not balancing the 
loss of energy from the radiating process, is not enough, 
and cannot be enough, regardless of any other feature 
(amplitude, frequency, polarization, etc.), to give rise to 
a strict stability of the orbits around the nucleus. Cer- 
tainly, as it already had somehow been made apparent 
after the work of Puthoff [2l| , the demand of strict com- 
pensation of radiated and absorbed power leads naturally 
to quantization (of AM, but not, at least simultaneously, 
on the values of the orbital radius and angular veloc- 
ity) as a necessary condition for stability. In no way, 
though, it constitutes a sufficient one, and the presence 
of a mechanism such as the one we will show here will 
turn out to be a necessary condition for that "sufficient" 
implication. Insisting on this issue a bit more, Puthoff's 
impressive derivation of the Bohr's ground state energy 
assumes equilibrium with the ZPF random background 
as a main hypothesis. He does not provide a real mech- 
anism through which to establish this equilibrium. 



C. What may be missing 

By radiative loss, a wider orbit will gradually converge 
to the "stable" one, and it is also true that, on aver- 
age, a narrower orbit will gradually grow, by absorption 
from the background, converging back again to that priv- 
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ileged one. But the key word here is "on average" . This 
average is equivalent to considering that the system is in 
an "statistical equilibrium" . Nothing, nevertheless, says 
how this equilibrium is achieved, or why instantaneous 
mismatches in radiated / absorbed power would yield "the 
right reaction" in the system. For instance, the excess of 
absorbed power in a "suboptimal" orbit could well make 
the particle return to the optimal one, or could well, re- 
calling that the ZPF is completely stochastic, make it 
collapse even faster towards the nucleus. 

From the point of view of a system theoretician, the 
extra element still needed in this picture is quite clear: 
we need some mechanism that allows the system to ac- 
cept the energy of a perturbation, store it, and then give 
it back again in the form of emitted radiation. This per- 
turbation would be associated with an instantaneous mis- 
match between the emitted (radiative loss) and absorbed 
(background) radiation. This kind of mechanism could 
now potentially allow an electron to orbit for an indefinite 
period of time around the nucleus, on a particular circu- 
lar, elliptical or more complicated trajectory that does 
not degrade over time. Naturally, at least in an average 
over a certain scale of time, radiative losses and absorp- 
tion from the background must compensate, if the orbit 
is to be stable. We believe that this kind of possibility 
is not at all generally acknowledged within the physics 
community. 

Aside from all that, the allowance for the system to dis- 
sipate is a strictly necessary condition for the appearance 
of what we call "attractors". This appears very clearly 
from energy conservation considerations. Indeed, the re- 
lation between radiated and absorbed power constitutes 
itself another feedback loop (FL), this time on average 
values. This is why the new dynamical FL that we will 
present here is so relevant: it applies to instantaneous 
values of the dynamical variables. 

Very graphically, Cole and others' previous work il- 
lustrates a background for these ideas. For an electron 
moving in a near circular orbit, under the effect of a cir- 
cularly polarized plane wave, normal to the plane of the 
orbit, we find (26j : "The result is a constant spiraling in 
and out motion of the electron, with the spirals grow- 
ing larger and larger in amplitude, until finally a critical 
point is reached and the decay of the orbit occurs" . For 
elliptical orbits, "this behavior also occurs for more gen- 
eral, but more complicated, elliptical orbits, where now 
an infinite set of plane waves is required to achieve the 
same effect, where the plane waves are harmonics of the 
periods of the orbit" . We need some kind of mechanism 
to force the electron to hold to its initial orbit. 

How can we do such a thing? We must find an inter- 
nal degree of freedom, which is the one that will store 
the energy. This degree of freedom will need to be able 
to "communicate" , to exchange energy with the "other" 
degree of freedom, the orbit, the one that is subject to 
external perturbation. The particular way in which this 
would result in a stabilization mechanism of the value of 
both is more a matter of the theory of dynamical systems. 



We will finish this little dissertation saying that spin, a 
"classical" spin, the rotational AM of the particle, seems, 
in a first guess, a perfect candidate for that "internal" de- 
gree of freedom we seek. Nevertheless, the interpretation 
of either the electron or the nucleus as a rotating uniform 
charge distribution leads to some difficulties with special 
relativity that were first encountered some eighty years 
ago, and will force us to be more concrete in our choice: 
we propose an interaction between the orbital degree of 
freedom and what we have called "residual spin" (RS, the 
amount of AM coming from the orbits of quarks inside 
the proton - if we had more than one nucleon, we could 
associate this RS with orbital AM of the nucleons inside 
the nucleus). 

As we had already hinted before, the picture is not 
complete, nevertheless, until we include some kind of 
background, because the system must be given a source 
to be able to compensate the radiative losses (dissipa- 
tion). Otherwise there is no way of making the system 
stable. But, still, how do we connect, in the classical 
formulation, those two degrees of freedom? 



D. Electromagnetic interactions within the 
hydrogen atom 

The quantum atomic Hamiltonian for the hydrogen 
atom consists of several terms. If we establish a relation 
with its classical counterpart, it is clear that the first two 
terms represent the coupling of the charge of the elec- 
tron to the scalar (electrostatic term, electric field) and 
vectorial (Lorentz term) part of the potential. These two 
terms configure the "coarse grain" atomic structure. The 
electron and proton lack on dipolar electric moment, the 
following (in order) non vanishing moment is their dipo- 
lar magnetic moment, but straight interaction (dipoling) 
of these does not appear as a primary term in the quan- 
tum development (the term of the kind a ■ a is not very 
relevant). 

Another term finally appears in the quantum formu- 
lation, giving rise to the "fine" structure of the atomic 
levels, namely the one known as LS coupling. At least 
two routes to this term are well known. The first of them 
departs from ordinary non-relativistic QM, and reasons 
semi-classically in the following way: the electron "sees" 
the proton orbiting around, and therefore "feels" a cer- 
tain magnetic field created by this current. Different ori- 
entations of the spin of the electron in this field yield the 
LS contribution, an additional factor 2 being introduced 
to account for the so-called Thomas precession. The sec- 
ond route to the LS term results from the combination 
of the Dirac equation with the principle of minimal cou- 
pling substitution. It is, therefore, fully relativistic, as so 
is the Dirac equation. 

All those terms are related either with coupling of fields 
to a point-like charge, or coupling of a field (a dipole field) 
with a dipole. All these interactions are somehow more 
"primordial" than the one we want to uncover. They 
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leave untouched the magnitude of (magnetic) dipoles, its 
action therefore reduced to a change of their orientation 
is space. 

E. Interaction between two systems of charges 

Classically, two magnetic dipoles of equal sign tend to 
anti-align. The magnetic field created by one of them 
causes a torque on the other. Moreover, the electric 
charge associated with the dipole also gives rise to a 
force, through the Lorentz law, F = q~E + v A B, with 
its purely electrostatic and magnetic terms, respectively. 
These two effects completely determine the interaction 
between two charge distributions, at least whenever mo- 
ments of higher order are all vanishing. Specifically, the 
first one of those effects, the classical dipoling interaction, 
would introduce an spin-spin term in the corresponding 
quantum Hamiltonian, where the second would stand for 
an "orbital" term, independent of the "spin" (in this case 
the classical self-rotation of the particle experiencing the 
field created by the source). 

In any case, our picture is more general to that of 
two charges with magnetic dipolar moment interacting 
electro-statically, through the Lorentz force, and also 
through the dipoling interaction. An "extended" dis- 
tribution of charges has an associated magnetic dipolar 
moment, but it is also capable of feeling the "gradient" 
of the field. This gives rise to some other effects that, 
had we considered an ideal dipole, would have remained 
unapparent. These new effects, as we will later show, do 
arise exclusively from the assignation of a spatial exten- 
sion to the charge distribution, that should for simplicity 
be regarded as a solid rigid (its movement being thus re- 
ducible, for any instant of time, to a pure translation with 
the velocity of the center of mass and a rotation, in this 
case around an axis of symmetry). In this last situation, 
the calculation of forces and torques no longer depend on 
the value of the fields at one point (the center of mass) 
but they involve the evaluation of an integral, where the 
action over each differential "element of charge" must be 
computed and added. 



F. A previous stage 

In Sec. lIIII initially, we focus on a somewhat simplified 
situation, and through a graphical argument, we expect 
to convince the reader that, under the action of a dipo- 
lar magnetic source, a (spherical) distribution of charge 
with non-vanishing second order moment (a spatial "ex- 
tension") experiences the following: 

(i) a net torque applied on its center, as a result of 
its "orbital" (center of mass, collective) velocity, in the 
plane normal to the dipoling source, 

(ii) a net resultant force, as a result of a non-vanishing 
AM of self-rotation, in the direction that joins the center 
of mass with the the dipoling source. 



In Sec. lIVI and subsequent ones, we do an already quan- 
titative analysis of the 2-dimensional problem. This sit- 
uation corresponds to what we have called "planar" ini- 
tial conditions: both the source magnetic dipole and the 
AM of the distribution are initially aligned, and remain 
like that because the initial velocity is strictly normal to 
them. Our intention is to be completely systematic, and 
to give quantitative expressions (although symbolic) for 
the two effects already introduced in Sec. IIIII 

Later, in Sec. IIV J[ we will complete our analysis al- 
lowing for general initial conditions. We will find out 
that the new contributions are (as expected) perpendic- 
ular to the ones in Sec. IIVI This is a highly satisfactory 
result because it means that feedback loop (FL) that we 
already had foreseen in IIVI is not destroyed by the new 
terms arising from the new freedom in the initial condi- 
tions. Sec. IIXI stands for an interpretation of the former 
results, from the point of view of the theory of systems 
dynamics (it is perhaps here convenient to say that the 
author has some background in this field). The theory 
of systems dynamics provides us with tools to determine 
necessary and sufficient conditions for the existence of 
stable "orbitals" in the classical configuration space of 
the system, without the need to address more specific 
calculations that in principle would lay beyond the scope 
of this paper. Something that may be clarifying is that 
we would be using here the term "orbital" in a much 
broader sense than just a kind of stationary /stable orbit: 
we mean just a portion of the configuration space of the 
system, where, once it is placed inside, there is a very 
little probability of leaving it. 



G. Model for a realistic hydrogen atom 

In Sec. |Vj we face now a slightly modified scenario: 
now the inertial system is anchored to the center of the 
charge distribution, and the point-like source is moving. 
This choice is, as we say, necessary to make the analogy 
with a real atom possible, where an almost massless point 
particle orbits around an almost static (because of its 
almost infinite relative mass) distribution of charge, the 
nucleus. 



H. A feedback loop (FL) 

Two variables "A" and "B" are under a FL if, upon a 
sudden change on "A", this influences "B" changing its 
value, and finally the change in "B" modifies again "A" . 
For instance, a time- varying electric field and its associ- 
ated magnetic field evolve under the action of a loop of 
influence. This kind of thing becomes apparent when we 
decouple the Maxwell equations for a perturbative cal- 
culation. Another example is the relation between the 
current of a coil and the charge of a capacitor in an os- 
cillating LC circuit. In our case, the "loop of influence" 
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in which we are interested relates the following variables: 

v<°' ^>t^u;^F^v( \ (1) 

where is the "orbital" velocity of the center of mass 
of the charge distribution, r is a torque applied at its 
center of mass, u) is an angular velocity vector, expressing 
a rotation around an axis of symmetry. By F we mean 
a resultant force, applied on that same center of mass. 
There are some subtleties regarding the definition of w, 
that we will treat later; for the moment we are content 
just to provide the reader with some "flavor" of what we 
are dealing with. 

Perhaps it is already convenient, even at this initial 
point, to enter a bit more into detail: actually, our equa- 
tions for the (2-dimensional) comprise two simultaneous 
FLs. The first of them will be responsible for the station- 
arity of the main orbital movement, and, once both the 
radiation losses and the interaction with a background 
are included, for its stability. Meanwhile, a second FL 
will introduce a secondary oscillation (2nd-Osc) that will 
allow us to provide a feasible explanation for excited 
states, hence for a full discrete and infinite, spectrum 
of (completely classical) stable states for the system. Fi- 
nally, when doing the extension of the analysis to 3D, a 
third FL appears, this last corresponding to the classical 
counterpart of the well known quantum LS. 

On the other hand, the association of these loops with 
particular "natural" frequencies of oscillation (which is 
crucial here) is proven (approximately, as the system 
is non- linear), via a frequency analysis of the linearized 
equations of the system. 

I. Linearization 

Given the dynamical equations of the system, what we 
do first if to look for stationary trajectories: if a suitable 
choice of dynamical variables is made, we can then as- 
sociate them with (classical) eigenstates of the system. 
We know, from system theory, that any harmonic func- 
tion is indeed an eigenstate of the linearized system, and 
this will be of great use later (though we use, however, 
more physical arguments to derive our stationary solu- 
tions, this idea will be of great use). A second (and last, 
for now) step is to identify which of those stationary 
trajectories are stable. To prove stability, we linearize 
the system around a "stationary" point x st in the form 
x = A-(x— x st ), where A £ M. nxn andx € TZ n is a vector 
of dynamical variables. Now, at least for that linearized 
description of the system, a necessary and sufficient con- 
dition for stability is simply that the real parts of the 
eigenvalues of the matrix A must be negative. 

III. AN ELECTROMECHANICAL GAME 

We work on a purely electromechanical model. This 
model has little relation to a real hydrogen atom: we 



are going to calculate the action of the field created by 
a point-like magnetic dipole on a distribution of charge 
with zero (monopolar) and second (quadrupole?) order 
non vanishing moments. The source is attached to an 
inertial system. The charge distribution can move as a 
rigid body. A reinterpretation the results that arise from 
this previous treatment to a more physical situation will 
be the main task of subsequent sections. 

We present now a very simple, qualitative argument, 
that will surely shed necessary light over all that we 
will calculate afterwards. Let us consider a point-like 
magnetic source with non-vanishing magnetic dipolar 
momentum, hence an ideal magnetic dipole. This 
source is attached (hence it does not move) to the 
origin of an inertial reference frame that we will call 
RF = {x ,i/o, z }- This magnetic dipole is aligned 
with the z-axis of RF n , and hence, far away from it, 
creates a magnetic field B in the direction — Zo that will 
act on any charge that may be moving with a certain 
velocity with respect to RFq, but it is not affected 
by them (this is equivalent to consider, for instance, 
that its mass is infinite). Now, we also consider a 
spherical distribution of (positive) charge whose center 
of mass is instantaneously moving with velocity v in the 
plane OXY with respect to RF 0} and rotating around 
its symmetry axis in the direction zq. We invite the 
reader now to a qualitative examination of two different 
situations. We have singled out points A and B to make 
the argument clear. 

1) Figure 1. The distribution moves as a whole with 
a linear "orbital" velocity v^°^ with respect to the 
source. There is no rotation (oj = 0). As a result of the 

radial component of this velocity (the component 
in xi), forces on A and B (and also at the rest of the 
points of the distribution) are induced. Clearly, because 
\Ba\ > |Bb|, I Fx I > |Fb| and there is a net torque r 
applied at the center of the charge distribution. The 
tangential component of this orbital velocity, v^°\ also 
gives rise to a net force, but we do not represent it for 
clarity. 

2) Figure 2. The distribution has no orbital motion now, 
but rotates around its center of symmetry. As a result 
of an angular velocity oj ^ (measured with respect to 
the inertial system RF ), the Lorentz law induces forces 
in points A and B, and all the points of the distribution. 
Again, because |B^| > |Bb|, \Fa\ > \Fb\, and (if we 
integrate for all points), there is a net resultant force 
acting on the mass center. 

Attention to sign consistency: from the graphical argu- 
ment, we see that a positive (counterclockwise) angular 
velocity to > yields a negative force on the radial direc- 
tion (F • x\ < 0, and a positive (outwards) radial velocity 

vi°^ > yields a positive (favors counter-clockwise rota- 
tion) torque (r • zq > 0) . There might be inconsistencies 
of sign in the rest of the document, but these two rela- 
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FIG. 1: Under the action of a magnetic source in the origin, 
an orbital motion (in the radial direction, vi"* 1 ^ 0) generates 
a net torque r =fc 0, due to the fact that |F^| > |F S |. 




FIG. 2: Under the action of a magnetic source in the origin, 
a self-rotational momentum oj ^ originates a net resultant 
force F_r / 0, modifying the orbit, due again to the fact that 
|F A | > |F B |. 

Hons are the base of the argument (we have checked their 
correctness in the dynamical equations). 

IV. QUANTITATIVE ANALYSIS IN 2D 

A. Reference frames 

As before, we consider an inertial reference frame RFo. 
The magnetic source stays attached to it, and the distri- 
bution moves. For convenience, we also consider a second 
system of reference RF\ , with its origin attached to the 
center of the distribution. The axes of RF\ are not rigidly 
attached to the distribution, so there can be a net angu- 



lar velocity lj for RF\ in respect to RFq. Indeed, we 
will let it rotate around an axis of symmetry of the dis- 
tribution. Moreover, xi, the x-axis of this system, RF\, 
will be always aligned along the positive direction of the 
vector that joins the point-like source and the center of 
mass of the distribution. We therefore define, formally 

Rorb = ~Rorb x l , (2) 

The auxiliary system RF\ will rotate instantaneously 
with angular velocity w with respect to RFq . Obviously, 
the fact that RF\ can rotate in relation to RFq clearly ex- 
cludes its inertiality. Nevertheless, we will express almost 
any vector in this frame RF\ . This is done for reasons of 
convenience. 

Very important: this vector system RF\ is only used 
as a reference in the mathematical sense, it is only a 
mathematical device. This means it is never used in the 
physical way: no force law (for instance, the Lorentz law), 
no mechanics (for instance, Newton's second law) is ever 
evaluated in this frame. What we can do, however, is 
express any vector (any position, velocity, AM, force or 
torque) as referred to its basis vectors Xi,yi,zi, which 
in turn are themselves functions of x ,yo,z . This use 
of RF\ will be very convenient to make the necessary 
changes that will allow to apply our picture to the more 
realistic situation of the hydrogen atom. 



B. Charge distribution 

For the present work, we will consider a charge den- 
sity p(t\) such that p(r{) > only for |ri| < R c , being 
the radius of the (spherical) distribution. We define a 
differential element of charge dq — p(ri) dv, with dv the 
differential element of volume. 



C. Purely kinematic considerations 

For each differential element of the charge distribution, 
we have a velocity that we will divide in two components: 

v = V W + v< fl »(ri), (3) 

where the superscript V will stand for "orbital" and 'R' 
will stand for "rotational" . For convenience, we will write 
all expressions in terms of the versors of RF\ . This poses 
no problem as long as we regard xi,yi,zi as functions 
of xo,yo,zo, and keep in mind RFo is our (only) inertial 
system. We will also include the dependence in ri when 
this is indeed present, as we do in v( fl '(ri). For the "or- 
bital" component we will barely write v^°\ as we will see 
it is not dependent in ri. We define, for each differential 
element of charge in the distribution: 

v<°> - Xl + « t (o) yi , (4) 
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with the components 



(°) (o) a ri 

v} = v\ ; xi = W( ) A R orb , 



dt ' 



(5) 
(G) 



where u>( ) is defined by ([5]) as the angular velocity of 
RFi with respect to RFq and 



o; (i j) A ri, |n| < i? c , 



(7) 



with ujir\ = W(r) z i (remember that in this section we 
restrict the initial conditions to the plane OXY), defined 
as 



J {R) 



•>(o), 



(8) 



and u! is the angular velocity of the distribution with re- 
spect to RF n (see [34])- Vectorially, U(r) = u) — W( )> 
with the three vectors directed in the Zi = Zo axis. It 
is important to remark here that u/ja, as defined above, 
is not an angular velocity as measured in RFq, in con- 
trast with w and W( Q ). This definition of uirm is justified 
for convenience, as it makes all our reasonings far more 
apparent. For clarification, and also because it will be 
useful, we add the following expression using spherical 
coordinates in RF\: 



J(R)ri sin 6»i [ sin^ixi -cos^iyi ] , (9) 



where r\ = and is positive for a counter- 

clockwise rotation. 



D. Fields 

We start by recalling that the electric field created by 
a point-like particle with charge —q (electric monopole): 



E(r ) 



4?reo |r | 3 ' 

but as a rule, we said we will express all in RF±: 
q ( r l — Rorb x i) 



E(n) 



47Te ki +Rorb Xi| 



(10) 



(11) 



with ro = i"i + R or b xi, with R or b the orbital radius. 
Acting on each differential element of charge there is a 
purely electrostatic force: 



and obviously 



dF e = dg(n) E(n), 



dF r , 



(12) 



(13) 



In the following we will face an analogous calculation with 
the Lorentz force. Here we have not developed (|12l) - (fT3"|) 
in a multipole expansion, as this would yield no special 



distinction among terms of a different nature. This will 
be, however, our main tool in the following calculations. 

Now, the magnetic field created by a point-like mag- 
netic dipole fjj = /xzi, if the source is sufficiently far to 
consider the observed field only has a z component: 



B(n) = -Bin) Zl , 

1 1 



I"l +Rorb Xi| 3 ' 



(14) 
(15) 



defining the quantity ft — ^fi, with /i > the value of 
the magnetic moment of the point-like source. 



E. Multipole expansion 

It is time to recall the first terms of a multipole expan- 
sion (around ro = Rorb, or equivalently, ri = 0\): 

1 



1 

|r | |ri - ( - R rb xi) 
1 



|ri + Rorb Xi| 

1 . {-Rorb xi) • ri 



Rorb Rorb 

1 Xi • r x 



Rorb R 2 or b 



1 



(16) 



|r X - ( - Rorb Xi)| 2 

1 

kl + Rorb Xl| 2 

1 2(-R orb xi) • ri 



/? 4 



+ 



1 2xi • ri 



R 2 R 3 

^orb ^orb 



1 



(17) 



ro 



|ri - ( - Rorb xi)| 3 

1 



ri + Rorb Xi| 3 

1 3(-R rb xi) • ri 



i? 5 



+ 



1 3xi • ri 



R: 



Ri 



(18) 



that we will now apply to the scalar B(ri), where B = 
B(ri)zi is the magnetic field created by the point-like 
source: 



B(n) = A 



1 3 xi • ri 



R 3 

^orb 



Ri 



(19) 
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On the electrostatic force 



If we now isolate the leading order, defining: 



We already settled the expression for the purely elec- 
trostatic force in (|12[) and (TT3"1) . That formula already 
accounts for the inclusion of higher order moments for 
the charge distribution. None of these terms, aside from 
introducing corrections to the overall magnitude of the 
resultant force F e , adds any other effect from our point of 
view. In fact, there is no component of F e in the tangen- 
tial (yi) direction, neither there is any kind of resultant 
torque over the center of mass of the distribution. For 
these reasons, we will not hereafter refer to this electro- 
static term F e , if it is not strictly necessary. 



G. Magnetic forces 

From here on, we concentrate on the forces arising from 
the "second term" within the Lorentz force, that of the 
kind qv A B. Acting on each differential element of 
charge, we have a force: 



dF r , 



dg(ri) v(ri) A -B(r x ) z 1 



dg(ri) 



,(«) 



a -B(n) Zl . 



(20) 

where we have applied our decomposition in ([3]) , and now 
we can write: 



dFW + dF<*>. 



(21) 



defining: 



dF$ = dq(rt) vWa-%) z 1; (22) 
dF™ = dg(r x ) v< fl >(ri) A -B(r x ) z v (23) 

1. Orbital component 

We start by analyzing the first of those two terms (a 
force from an "orbital" origin): 

dF$ = rfg(n) V W A-S(n) z x 

= dg(n) v< 8 ' A 

1 3x! • r! 



B 3 



i? 4 

^orb 



zi, (24) 



and now, integrating for the whole distribution: 

F<°> = | dpW, (25) 

and taking into account ([T91 , we clearly see the first term 
(electric monopole) already gives the leading contribution 
to the integral, and therefore: 



F# = J dg(n) v<°> 



(v (o) A zi) 



A /t ■ 



-"-orb 



Zl 



7? 3 



dg(ri). 



(26) 



(v( o) A zi 



Qo = J dq{v{), 



(27) 
(28) 



where Qoj & zero order moment of charge around 0\ 1 i.e., 
Qo = Q, with g > the total charge of the distribution. 
Let us now use Q and give a more detailed expression: 



p(°) 



Ft 3 



(29) 



though, in this case we could simply use the full order, 

PW = |B(r = i? o ,fc)| 
- |S(ri=0)| 



(°) (o) 

^* x i - «r yi 



Vt (o) xi-^) yi 



(30) 



2. Rotational ( "spinning") component 
Recalling (O, we set now the focus on v( R \ 

dFW = dg(ri) v< fl) A -B(n) z x 

= dg(ri) (w( R )Ari) A -B(ri) zi. 

We can calculate a bit more: 

dF^ = dg(ri) zi A ri A --B(ri) zi 

= dg(r x ) W( fl )B(ri) z x A zi A r[ h) 

= dg(n) u^BW (-r< h) ) 

= - dg(ri) w (fl) B(ri) r£ \ 



(31) 



(32) 



where each time we do zi A on the left we "rotate" a 
2D vector by an angle ir/2 around the z. We have also 
introduced, for a horizontal projection of a vector, the 
following notation, for any vector A: 

AW = (A-xi) X! + (A-yi) yi, (33) 
A<*> = (A-zi) Zl . (34) 



We return to the integration now 

p(R) - f .vpCR) 



(35) 



we can see the first term (1191) gives rise to a contribu- 
tion of first order to F^f (be aware of the last factor 

ix !), thus vanishing (as the dipolar electric moment of 
the distribution does). The second term in (|19[) gives rise, 
however, to a non vanishing contribution (as the second 
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order moment does not vanish), due again to that last 



factor r[ h \ We can then write: 
- J dg(ri) u w fi 
dq(ri) uj (R) ji 

= / dq{vi) 



3xi ■ n 



(~r[ h) ) 



3xi ■ n 

RL. 



(h)|2 



Xl 



3 r 



(M|2 



R orb 

3\v[ h) \ 2 



Xl, 



(36) 



where we have used R or b = Rorb xi. We now truncate 
to leading order, and solve the integral: 



U orb 



(37) 



where Q2 is a second-order axial momentum around the 
Zi-axis: 



Q-2 



dq(n) \r\ 



3. Remarks 



(M|2 



(38) 



We remark that both d-Fm ^ , dFffi , under these 2D ini- 
tial conditions, only have components in the plane OXY: 

F m = J ( dFM + dFCR) ) = | dFW, (39) 

and we have also seen the first non vanishing term in both 
dFm \ in (|25|) . and dFm , (1351) . belongs, respectively, to 
zero and second order (desarrollar esto!!!). 

H. Magnetic torque 

We analyze now the torque, that we previously intro- 
duced, from a very simple but intuitive graphical argu- 
ment (reference to graphics). As a summary of this sec- 
tion, let us say we will show that this torque does indeed 
comes exclusively from the orbital movement, as we had 
already foreseen through the graphical argument. We 
define: 



.(°) 



-(R) 



dr$, dr%=T X A d¥t\ (40) 



JdrW, dr(f)=n A dF£>, (41) 



and obviously t t 



» 



1. Orbital component 
Using (f22j) we have 

drtf = r i A dF% 

= ri A dg(ri) [ v (o) A -S(ri) zi . (42) 



Now, with v^) = !)f xi + because our planar 



initial condition, we do 
dr$ = ri A dg(n) 



xi + v[ o) yi 



A -B(ri) zi 



dq{vx) Bin) n A (-«l°J yi + v[ o) xi), 



(43) 



Also, due to reflection symmetry around the OXY plane, 
all contributions to t$ in that plane vanish when inte- 
grated: 



dr# = | (drW) (z) , 

and using (|4"2"j) we do 

dr# = ri A dg(n) [ -v<°' A -B(n) z x 
= v (o) A dg(n) ri A -S(ri) z x 



(44) 



= -v^ A dg(ri) A B(n) z l5 (45) 

and applying the multipole expansion for the magnetic 
field 



dr^ = - V W A dg(n) 



( A p, 



1 3xi • ri 



Ri 



Ri 



(46) 



With the "additional" factor , the first term of the ex- 
pansion goes from zero to second order. Clearly it gives 
no net contribution to the integral. The second term goes 
from first to second order, and it is the leading contribu- 
tion. Keeping this leading term and truncating: 



dr$ « -vW A dq{rx) r\ h) A /t 



v<°> A dq(n) r\ n > A A 



3xi ■ ri 

i? 4 

-"-orb 

3xi ■ n 

p4 



Zl 



= vW 



v^A fidqin) 



3\r[ h) \ 2 cos^! 
/? 4 

x ( sin?/>iXi — cos^iyi ) 



Q|_( h )|2 
-"orb 



x ( sin ^1 cos ^1X1 — cos 2 V'iyi ), 

(47) 
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We could not avoid having to invoke explicitly some co- 
ordinates here, (f>i being the azimuthal angle in system 1. 
The term in sin <pi cos 4>i clearly vanishes if we integrate 
for ip € [0, 2ir], but no the one in cos 2 ipi, and finally: 



-0) 



v<°> A A^(ri) 



A 



p4 
^■orb 



3|r^| 2 cos 2 ^ 

r>4 

»|2 



yi 



dg(n) |ri"^cos^i yi 



yi 



i? 4 

^orb 

^Q 2 a (o) 3 AQ 2 

—3 — yi A = 4 



(*) 



(48) 

in the ^-direction, as expected, and where 

Q« = J '^(nHrf^cos 2 ^, (49) 

is a second order axial moment of the distribution (we 
have to check this calculations). 



2. Rotational component 

Now we analyze dr\n^ , and we prove that TnP = 0, 
i.e., the only contribution to this torque of magnetic ori- 
gin comes from the orbital term. We see: 

dr£> = n A dFW 

= ri A dq(rx) ( v (fl) A -B(rt) z x ) 

= dq(ri) n A (w (R) A n) A --B(n) z x 

= dq(ri) n A W(ji)Zi A B(n) z 1 A n 

= dq{r x ) n A w (Jl) B(ri) A (-r^) 



dqfa) u W B{n) rf 3 A 

-dq(ri) w (fl) S(ri) } 
0. 



A r 



-1 
(ft) 



Therefore, 



0. 



(50) 
(51) 



It is worth to remark this last result: the torque on 
the z-direction (to be applied on the center of mass of 
the distribution of charge), comes exclusively from the 
"orbital" contribution, that contribution with origin in 
the "orbital" AM of the distribution as a point particle 
(therefore, represented by a center of mass). 



I. Summary of results for the 2-dimensional 
problem 

The following shows the first contribution in the elec- 
trical multipole expansion for each of the forces or 



torques. To second order (Q), no other forces or rota- 
tional momenta arise. 



Order 


pW 


■*- m 


r (o) 


T («) 

• m 


Zero (M) 


Yes 


No 


No 


No 


First (D) 


Yes 


No 


No* 


No 


Second (Q) 


Yes 


Yes 


Yes 


No 



A table of dependencies may be of use, too, 



Variable 


■If m 


■*- m 


_(«0 _(«) 

' m 1 m 


|-Ror6 1 


Yes 


Yes 


Yes 


U{R) 




Yes 




(o) 
Vr 


Yes 




Yes 


C°) 
Vt 


Yes 







J. 3D initial conditions 

We have a 3D extension of the problem. Here, xi 
still goes in the radial direction of the orbit, yi in the 
tangential direction and Zi is always the instantaneous 



axis for 



We have now: 



v (o) = ^ o) Xl 



(°) 

vl yi 



and 



>(R) 



(w(m)r xi + (w(m)t yi + [u(r)) z z 1; 



with, now, vectorially (recall Sec. HVCj) . 

U(B.) = W — W( ), 



(52) 

(53) 
(54) 



Needless to say, we can always decompose the dynamics 
of the system in this way. Now we complete the sys- 
tematic analysis of the dominant terms for the forces, 



F (°) F (R) 



m , m , and torques, r m , , m 



Tm . First we calculate the 



new contributions to the forces: 

A(dF&>) = dg(n) 4 o) Zl A -B(n) Zl = 0, 



(55) 



so there is no new contribution to F„ , and 

A(t£P^) = dg(n) («(«)) w A n A -B(n) Zl 



dg(ri) S(n) (w ( m) ( ' l) A n A zi 



A(dF^) 



(flh(z) 



(56) 



This is rewarding because the new contribution is normal 
to the one we already have. We can also write 



A(dP^) = dg(r 1 )B(r 1 )(u; (Ji) )W A v\ 



(h) 



(57) 
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where r^'*^ 2 ^ stands for a rotation of r[ h ^ by an angle 
ir/2 around the z-axis. We would have the integral 

A(F m ) = J A(dF rn ) 

= J A(dF^) + J A(dF<«)) 

= J A« fl ))W = A(F<«))W, (58) 

where we have just identified that F^f ' only has a compo- 
nent in the z-direction. Now we see the new components 
for the torques. Clearly from (|55[) we see there is nei- 
ther any new contribution to the part of the torque that 
originates from the "orbital" movement: 

A(dT^) = n A A(dF&>)=0, (59) 

and for the rotational part, using that A(dF^) = 
A(dFm we can write 

A(dr^) = n A A(dF^) 

= n A A(dF^) w 

= r^A A(dF^) (z) = (AdrW)W, 

(60) 

again with the rewarding result that this new contri- 
bution A(dTnP) is a "horizontal" vector (vector in the 
plane OXY) and hence normal to the contribution we 
already had for the "planar" initial conditions. 

Again we attempt to summarize everything in a table. 
These results are generalized to any -planar or not- initial 
conditions. First line says whether a net contribution 
exists. Second line establishes a dependence. Last line 
says which is the leading order: zero (M), first (D) or 
second (Q). 



* m -T m 


■t m 


-f m 


• m i m 


(R,h) 
' m 


Yes 


Yes 


Yes 


Yes 


Yes 


(o) 

Vr — 






- vl o) 


("(R)) W 


M 


Q Q 


Q 


Q 



On the other hand, the new terms should correspond 
to what is known as the classical counterpart of the quan- 
tum LS or spin-orbit coupling. This interaction is also a 
stabilizing one, as it tends to keep parallel the spinning 
axis of particle and the axis of the orbit (the direction 
of the dipoling source). Nevertheless, for the moment we 
will only pay limited attention to questions regarding the 
problem in 3D. 

V. TOWARDS A REALISTIC MODEL OF THE 
H-ATOM: RELOCATION OF MASS 

We had begun presenting an idealized electromechani- 
cal game. Now we establish a bridge from that situation 



to a realistic model for the classical H-atom. This implies 
a relocation of the inertial system, and a proof that the 
equations of the system remain the same in the new sit- 
uation. Later, the value of the parameters in the model 
must be adjusted to resemble the actual charges, masses, 
etc of the hydrogen atom, but this will not be, for con- 
venience, done yet. 



A. Relocation of mass 

In this section we do a relocation of the inertial frame 
of reference, RFq. To indicate this we introduce the no- 
tation RF \( n \ in contrast to RF \( old \ If our situation 
is to resemble reasonably the real hydrogen atom, a natu- 
ral choice is to attach RF$ | W to the center of the proton 
(the charge distribution) , taking advantage from its much 
higher mass in respect to the electron. Still, the proton 
can "rotate" around zq (on the planar problem), with 
angular velocity 

w| (n) = W n) + W(o) |("\ (61) 

with respect to RFq \ ( n > . This definition is consistent with 
the one already given in ([5]), so, analogously to what hap- 
pened there, u>\^ and W( )|^ both have an "inertial" 
meaning (they are both angular frequencies measured 
with respect to an inertial reference frame RFq), while 
W (-R)| ^ does not (but is nevertheless a very convenient 
dynamical variable from the point of view of the equa- 
tions and our whole argument). Of course, we are only 
dealing here with the 2-dimensional situation: for the 3D 
problem, the same would hold but vectorially this time 
(it is always possible to decompose the (instantaneous) 
movement this way): u>|w = + o;^ ^|( n ). 

Again, an auxiliary system i?Fi|( n ) is of good use. 
We choose to attach it, as we did before, to the cen- 
ter of the distribution, and also choose the same ori- 
entation, so that i?Fi| (n) = i?Fi| (oW) (hence we simply 
write RFi), but this time, also, the origin RFi\( n ' and 
RFq\^ is common. On the other hand, the point-like 




FIG. 3: Location of systems RFq and RF\ in situation (A) 
and (B). In (A), origin of RFo coincides with the position of 
the magnetic point-like dipole /x, while in (B) both origins 
coincide. 
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magnetic source (an electron, here) will orbit with veloc- 
ity v (°)|(") = v t (o) |(") + vf o) |(' l) around the proton, and, 
again in this situation, the components can be expressed 
using the versors of RFi, in the way: 



,(<0 



l(n) = »r |(n) x l, 
v t l(n) - v t l(n) Yl- 



B. Formal equivalence of the equations 



(62) 
(63) 



The problem of a moving magnetic dipole is a classical 
one in electromagnetism. Due to Lorentz covariance, a 
magnetic dipole fi that moves with velocity v in respect 
to an inertial frame is seen, by an observer in that inertial 
frame, as an electrical dipole of value p proportional to 
fi A v. For more details in this result, we can cite , for 
instance, [3l[ (18-5, page 334). This electric "effective" 
dipole creates, at the point i"o = R or b xi + rj., an electric 
field 



El 



(n) OC 



M A v l(n) 



\Rorb Xi + ri| 3 \R orb Xi + Ti| 3 

therefore exerting a force 

M A v| (n) 



dF\ 



(n) 



dq E| ( „) cx dq- 



Rorb Xi + ri| 



, (64) 



(65) 



on every element of charge dq of the distribution. Now, 
we recall that in our previous situation we had, suffi- 
ciently away from the source, 



B 



(old) 



Mo 



and therefore 

rfP| (oW ) - dq "-' R3 - = ,/,, 



47T |i? orb xi + ri | 3 ' 

v|( oi <i) A (i 



(66) 



i? orb xi + n 



(67) 



But now is the source that moves in respect to the ele- 
ment of charge dq. Clearly, v|r n ) = —v\( id)- Moreover, 
for clarity we also state: 



v (o) | ( „) = -v<°>| (oW) , 
vW| (n) = -v<*>| (oW) , 



and also clearly, 



(n) 



U! 



(n) 



(68) 
(69) 

(70) 



But it does suffice to perform the substitution v|( n ) = 
— v l(oid) in ([65|l . and so we have, again, 



, v|( oW) A/x 
dF (n) oc dg— ■ ^, 



(71) 



This is the same as we had in (I57|) and so it was exactly 
what we were looking for. 



Therefore, at least for the term that depends on a 
wedge product on the velocity vector, the expression 
of the force F on each element of charge is completely 
equivalent (modulus a certain constant) to the one in 
our previous situation, the first electromechanical model 
we presented here. Moreover, all previous expressions 
are applicable (modulus a possible constant), and no fur- 
ther changes needed, as we take advantage here of the 
fact that all of them were referred to RF\ , that has not 
changed in the new picture (using this trick has saved us 
a lot of calculations). 




1 






















— > 



FIG. 4: Direction of the force in situation (A): dF = dgvAB, 
and (B): dF — dqE. From RFq, the magnetic dipole (i is seen 
as an electric dipole p, source of the corresponding electric 
field. 



VI. DYNAMICAL EQUATIONS (2D) 

We have seen that the equations of the system are in- 
variant about whether it is the point-like negative charge 
or the extended distribution (of positive charge) that is 
moving, if we do the convenient redefinition of velocities. 
Using that fact, we present a first set of equations, that 
will be later enhanced by the inclusion of the radiative 
correction as well as the absorption from the background. 
This second step is done, anyway, only in an "approxi- 
mate" way that is nevertheless enough for our purposes. 



A. Dynamical equations in absence of rad/abs 

In the following equations, all quantities (q, E\, B±, 
etc.) are defined positive. In absence of dissipation and 
absorption and for the planar problem, the dynamical 
equations of the system are, to leading order, the follow- 
ing: 



fa (o) ) 2 . 

Rorb 

- qE x 



qBi v, 



7n m 

(o) 3/iQ 2 



• xi 



R 



4 

orb 



(72) 
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(o) 



4- v(R) 



■ yi 



qB\ v. 



M- 



2 W(R) 



-(°) . 



Zl 



zi - r (o) 



yi 



Zl 



3 A „(o) 

i? 4 r ' 

-"-orb 



(73) 



(74) 



with M 2 a mass second order moment or "inertia" mo- 
ment. In the first of the former equations (|72p . we take 
into account the fact that i2.Fi is rotating with respect 
to RFq, and that we have defined this radial velocity as 



,(°) 



,(°) 



? (R) 



component in yi. Also, we have to notice E\, B\ do de- 
pend on R or b (see below). To the former three equations, 
we must add the obvious 



i2o 



,(°) 



(75) 



and they must also be supplemented with the values of 
the field evaluated in ri = 0. With E(ri) = — E(ri) xi 
and B(ri) = — B(ri) zi, we have 



|E(ri)| = 
|B( ri )| = 



q 



4^0 R 2 orb 



MO Me 

4tt i2i 



(76) 
(77) 



where /i e is the magnetic moment of the electron. See 
PT|) - (ll5p . It is important to note that these fields enter 
the former equations with the values in the point ro = 
R or b (i.e., ri = 0). To leading order, then, any other 
information about the distribution is already contained 
in g, Q2 and Q^- The first moment Q\ — vanishes, 
as all charges in the distribution arc of equal sign, and 
moments of higher order are not considered. 



B. Inclusion of rad/abs 

How do we include the radiative correction? From the 
point of view of the equations, it would suffice to include 
a stochastic term in (|72l) . (f73|). and possibly in (|74)) . This 
stochastic term represents the difference of loss and ab- 
sorbed radiation in each instant of time, and its mean 
value is zero for a stationary orbit (an "attractor"). For 
non stationary orbits the mean value of this stochastic 
term would not be zero. The inclusion in (f?2"j) . (fT3"f is 
quite obvious: we must allow for an energy loss through 
a fall of the orbital radius and a decrement of the modu- 
lus of the tangential velocity. The term in ([74]) is not so 
obvious, and we will extend on this later. On the other 
hand, a more rigorous treatment is not necessary for our 
purposes, at least for now. 



1. Rad/abs for a circular orbit 

Puthoff calculated (P a b) (orbit) assuming statistical 
equilibrium with the orbital degrees of freedom (these are 
the two spatial coordinates in the plane, oscillating with 
frequency , so therefore they can be seen as two one- 
dimensional harmonic oscillators in quadrature) . Follow- 
ing 21]. we had 



(Prad)( 



„2 e>2 



circ) 



(78) 



which is directly obtainable from the Larmor formula 
with acceleration 



•(c) _ [Vt ) 



R, 



orb 



^(o)Rorb, 



and 



(Pab}( 



q 2 fu*)f 



circ) 



(.0) 



67T60WC 3 ' 



(79) 



(80) 



But here we deal with an instantaneous basis: the par- 
ticle will suffer the action of a field with an stochastic 
instantaneous value (the ZPF), and will also loose en- 
ergy whose (expectation?) (instantaneous) value is given 
by the radiation term (dependent on the instantaneous 
velocity). 



2. Rad/abs for the radial component 

To first approximation, we can introduce a correction 
just in (|72|) , through the inclusion of an stochastic term 
Xr, with an expectation value which depends on i2 or b, 
and must change sign around the "stable" value R or b- 
This i2^' b is the value of i2 or b for which the balance, on 
average over a cyclic orbit, of loss and absorption takes 
place, i.e., an equality between eqs. (fBTjf and ([75)) holds. 
So, 



(Xr(Rorb)} < 0, 
(Xr(Rorb)) > 0, 



Rorb > R-orl 
Rorb ^ Rnr 



b- 



for example, 

{Xr{Rorb)} 



kr (Rorb 



Rorb), k r > 0. 



(81) 
(82) 

(83) 



For clarification, we have to say that the instantaneous 
values of P a b and P ra d are stochastic variables whose dis- 
tributions depend on the instantaneous values of i2 or b, 
vi°\ etc., and therefore, strictly speaking x(Rorb) — 
f (Pab — Prad), and (|83"|) is only justified for mean val- 
ues (which is what we have done). 



3. Rad/abs for the tangential component 

Moreover, the gain/loss of energy affects both the ra- 
dial and tangential components of the velocity. For 
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this reason, we consider a second stochastic component 
Xt(Rorb), satisfying, this time 

( Xt K (o) ))<0, v^^v^U (84) 
( Xt K (o) ))>0, v\ o) <v[ o) \ st , (85) 

and therefore, close to the point of equilibrium \ s t, we 
can write 

(Xt(«t (0) )>«*t(« t W U-«t (0) ), k t >0. (86) 

4- Rad/abs due to "spinning" 

Our study of the dynamical equations lead us to con- 
clude that we need another dissipation/gain loss for the 
degree of freedom represented by uj. The rotational move- 
ment of the distribution around an axis is also subjected 
to loss and absorption from the background. In the way 
we have done before, we write 

(X«(w)> « k u (uj st - u), K>0. (87) 

Naturally, for a realistic hydrogen atom we would have 
uj\ st =0 (the proton is quasi-attached to an inertial sys- 
tem: it can oscillate with respect to it, but the mean 
value of this oscillation is zero). We will comment on 
this later. 



VII. FEEDBACK LOOPS 

A. A primary FL (lst-FL) 

As a preliminary approach, we are interested in a "pri- 
mary" feedback loop" (lst-FL), 

4°) r => uj => F => vi°\ (92) 

where r and F are a torque and force, as we already 
said, mediating between the (orbital) radial component 
of the instant velocity and the self rotation of the charge 
distribution. We can also interpret this in terms of the 
orbital radius, 

Rorb => r => uj => F R orb . (93) 

Now, if we look at the equations (|89p and (|9ip. we can 
identify a FL with an odd number of "minus" signs. We 
now explain what this means. Equation (|90[) will be left 
aside for the moment, for the sake of clarity. Whenever 
there exist these kinds of FLs, the negative sign of one 
of them is a necessary and sufficient condition for the 
existence of some kind of stability (this is a well known 
result from the theory of dynamical systems). 

B. A secondary FL (2nd-FL) 



5. Complete dynamical equations 

We now include all the former corrections in the dy- 
namical equations of the system. First, for the radial 
component: 

mi)), ' + m 



Rorb 

-qEi - qBi v y t - — ^- u) {R) + Xr, 

U orb 



(88) 



and reordering (that term on the left...), 



mv 



(o) ~ - qEl - qB 1 v[ o) - m ^-L 



Rorb 

-j^U {R) + X r, (89) 



Secondly, for the tangential component, 



mv { t o) k +qB l «W + X t, 



(90) 



and to conclude for now, for the "spinning" of the distri- 
bution, 

M 2 U ( R) « + 4 o) + Xu- (91) 

U orb 



There is a "secondary" feedback loop (2nd-FL) in the 
equations, 



,(°) 



,(°) 



(94) 



that would prevent the existence of stationary orbits 
where vi (t) = 0,Vt. However, we will in principle dis- 
regard this effect for simplicity (simply ignoring the cor- 
responding term in the equations). Later, we will see 
that it is precisely this secondary loop the reason why 
stationary /stable orbits corresponding to a quantum la- 
beling n = are not present in the real spectrum, as so 
it happens for n = 1 for orbitals with net AM (there are 
no real orbitals for n = 0, only s-orbitals for n = 1, only 
s,p orbitals for n = 2, etc.). 

The presence of a natural frequency of resonance is 
directly linked, and proven, by the fact that all poles of 
the system (under linearization) are purely imaginary. 
Indeed, the ( approximate ) linearization of the system 
around a stationary point ( either a circular or pendulum- 
like orbit) yields a set of purely imaginary poles (in the 
frequency domain). This corresponds to a harmonic os- 
cillatory behavior. 



VIII. A HYDROGEN MODEL IN 2D 

Up to here our equations were completely general, now, 
to particularize the equations to the hydrogen system, we 
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must make the substitutions (check the correspondence 
in each equation): 



rn 
M 

Mo 



Tfl e , 711, p , 
Me, 

„ _ Mo 

^e - 4?r Me, 



(95) 

(96) 
(97) 



with, naturally, m e , m p the masses of the electron and 
proton, respectively, and fi e the magnetic moment of the 
electron. On the other hand, for a realistic hydrogen 
atom, (lo) = 0, i.e., the proton, with an inertial frame 
attached to its origin, will not move "on average", but 
can rotate on small oscillations around the axis of the 
orbit (at least in the circular case). This point is for key 
importance and we will return to it several more times. 



IX. 



AN APPROACH TO STABILITY WITH 
NET AM 



An overview 



In the first of the following subsections, IIXB1 wc be- 
gin by presenting a qualitative reasoning, that may be 
enough to convince the reader but that it cannot be con- 
sidered yet a strict proof. Then, we do some calculations 
to show that a stationary point indeed exists for the sys- 
tem of equations that describe the system. This station- 
ary point, however, can be stable, unstable or critically 
stable. In the first case, the system will answer any per- 
turbation that drags it (not very) far from the point with 
a reaction that tries to attract it again to the stationary 
point. In the second, any perturbation (no matter how 
small) will launch the system on a trajectory that will 
gradually distance it from the initial state. The third 
case, critical stability, accounts for a non-limited storage 
of energy, in the form or perpetual oscillations. 

It is the aim of subsection IIX II to discriminate between 
those three possibilities, determining which one of them 
occurs. This is done through a very clear and well defined 
mathematical condition: the real part of all eigenvalues of 
a certain matrix, involved in the dynamical description 
of the system, must be strictly negative. Nevertheless, 
through our qualitative approach we already have the 
certainty that, at least on a certain point, the system is 
indeed strictly stable (the first of the three possibilities) . 

All this, for the moment, applies strictly only to the 
trivial case of circular orbits in the 2-dimensional situa- 
tion, although the extension to the 3D space, as well as 
to more complicated cases, like elliptical orbits, can be 
already, and somehow, be foreseen from the base that we 
are settling here. Besides, we will not do any estimates of 
energies yet here, that will be a matter of future sections. 



B. Qualitative approach 

1. The lst-FL makes stability possible 

Now, if we look at the equations and (|9"Tj) , we 
can identify a feedback loop, that we already treated in 
Sec. IVII Al Whenever there exist these (negative) feed- 
back loops, the negative sign of one of them is a necessary 
and sufficient condition for the existence of some kind of 
stability (this is a well known result from the theory of 
dynamical systems). 

Following (|9"Tj) . we see that an increase in Vr°\ the 
radial component of the "orbital velocity" , causes an in- 
crease in U(r) (and therefore also in oj = w/r) + wt \). 
But, through (|59"1) . and increase in lj^ gives as a conse- 
quence a decrease in vi°\ We remind the reader here 
that lo and u>r \ are "inertial" angular frequencies, the 
first of them referring to a self-rotation of the proton 
(the residual part of the inner AM coming from the fact 
that it is a charge distribution, most of that (inner) AM 
-therefore spin - being carried by quarks as point-like en- 
tities), and the second referring to the orbital movement 
of the electron around. Meanwhile, has no inertial 

meaning and is just a convenient dynamical variable to 
work with. 

This stability could be in principle strict, in the sense 
of a strictly stable point, or could also be a saddle point, 
or/and, given the nonlinear nature (see below) of the dy- 
namical equations, it could give rise to more complicated 
stable "orbits" in the configuration space of the system. 
We use "orbit" in a somewhat broader sense than in pre- 
vious sections, meaning this time any collection of points 
in the configuration space of the system where the prob- 
ability to find it does not change with time. Stability 
implies that the system, confronted with a small pertur- 
bation, always responds in a way that tries to compensate 
that perturbation. The "smallness" of the perturbation 
sets it clear that stability is, in principle, a local concept, 
valid in a certain "environment" within the configura- 
tion space of the system. As we already said in previous 
sections, this perturbation would be, for example, a sud- 
den "push" , due to a mismatch on the emitted/absorbed 
energy. 

A simple argument could be this: let us suppose the 
system is in a circular orbit with Vr =0. If, by ac- 
tion of this perturbation, the system is "pushed" afar 
to a higher orbit (there is a positive Vr°\ and the or- 
bital radius R or b increases), the system reacts increasing 
tjJim, thus making ixf negative. Oscillations may hap- 
pen, and, if the system is allowed to dissipate (in this 
case, by radiation prevailing over absorption), the initial 
configuration will be recovered after some time. This be- 
havior can be summarized by saying the excess of energy 
is temporarily stored in that inner degree of freedom that 
we have named "residual spin" (residual oscillations of 
the AM of the proton, coming from the "orbital" move- 
ment of quarks inside, in our relativistic interpretation). 
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From there, it will be released as radiation in subsequent 
instants of time. 

If the effect of the perturbation is to make the orbit 
loose energy, the process is the opposite one: as the ini- 
tial Vr°^ is negative, to reach a lower energy orbit (lower 
orbital radius), the system reacts making u)/ R \ decrease, 

and this makes vi°^ increase again and turn positive. 
This is equivalent to say that the system compensates 
the sudden loss of energy drawing it from the "residual" 
spin degree of freedom, to which it will later be returned, 
after the systems absorbs it from the background. Of 
course, a sufficient quantity of energy should be dispos- 
able in the immediately following instants of time or this 
scheme could not work. As we have said, our bath of 
stochastic radiation will be the provider of that energy 
in subsequent times. 



2. Consequences of the 2nd-FL 

The secondary feedback loop v[°^ =>• vi°^ => v[°^ will 
prevent stationary orbits where vi°\t) — 0,Vt. How- 
ever, we will in principle disregard this effect for simplic- 
ity (simply ignoring the corresponding term in the equa- 
tions). Later, we will see that it is precisely this 2nd-FL 
the reason why stationary /stable orbits corresponding to 
a quantum labeling n — are not present in the real 
spectrum, as so it happens for n = 1 for orbitals with 
net angular momentum (there is no s neither p-orbital 
for n — 0, and no p for n = 1 either). 

3. Absence of rad/abs: a continuum of stationary orbits 

By a stationary trajectory we mean a closed, cyclic 
one, such that if the system is in a particular point of 
the it, it will describe on and on that same one, coming 
back to exactly the same initial condition after each cy- 
cle. We will also use the term "stationary orbits" . Later 
we will introduce the concept of "stationary set" as a set 
of points or initial conditions given that, if the system is 
initially in one of them, it will evolve without leaving that 
set. This generalizes the concept of stationary trajectory. 
Stationary trajectories or orbits can be stable or unsta- 
ble, depending how they respond to small perturbations, 
once the stochastic terms are included in the equations 
of the system. In absence of radiation, a perturbation 
causes a transition between two of our continuum of sta- 
tionary trajectories. 

4- Forcing rad/abs balance 

The following step in our logical development is the 
following: when imposing the rad/abs balance, the result 
is that a particular one of the former family of strictly 
circular trajectories is now singled out. It behaves as 



an attractor in the configuration space of the system. If 
we simulate the effect of the random background on a 
particle traveling in that privileged orbit, we will get a 
probability distribution extended to the whole space, in 
other words, an orbital. 

C. A set of stationary trajectories 

As we have said already, in this section we will disre- 
gard the presence of that 2nd-FL, i.e., we simply ignore 
the term connecting vf^ andvf in eqs. t8fy) or \89j) . We 
will show that there is a continuum family of stationary 
trajectories in respect to the dynamical equations of the 
system. The stationary points in this collection are not 
stable in the strict sense, though: a (small) perturbation 
will cause a transition, mediated by certain oscillations, 
between two of these trajectories. Later, the inclusion of 
the rad/abs balance will single out one of these trajecto- 
ries, also providing it with the feature of stability. 

With no need to get into mathematical analysis, we 
already now there is stability in the case when vi°^ = 
(hence, v { r o) = 0). If v ( r o) = 0, from ([73]), v ( t ° ] = 
and = too. This means that stationary circular 
orbits are given, in this circular case, by the equation that 
guarantees vi°^ = (again ignoring the term connecting 
vi o) and v[ o) in eqs. ([88 ]) or (|89 ]l . We have, from, ([72]): 

ip , n (o) , 3Ae<22 K (o) ) 2 (n8 , 
qEi + qBi v\ + — uj (r) = m e — , (98) 

where we have used m e , the mass of the electron, at the 
right hand side (it is the electron that revolves around 
the proton), where jl e — and where Q2 is, as usual, 
a second order charge momentum (expressing the fact 
that the proton is seen as an extended distribution of 
charges). We can recognize, on the left hand side, the 
three leading contributions on our multipole expansion 
of the problem: (i) electrostatic attraction, (ii) Lorentz 
force on a point-like charge, and the third one, (iii) a 
residual Lorentz effect on an "extended" particle. 

Now, with v[ qS> = bJt \ R or i, for a circular orbit, we can 
rewrite the former equation as 

qE\ + qBi R or b H —4 — W(_r) = m e R or b w ( )> 

"■orb 

(99) 

with the values of the fields depending on the radial co- 
ordinate Rorb, 

Ei = |E(n)|=+ q , (100) 
47re o R orb 

B 1 = |B( ri )|=+^- = ^-, (101) 

"•orb "'orb 

again with fi e = ^-fi e (fi e the magnetic moment of the 
electron), and where q > is the unit charge. We finally 
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arrive to an equation 



<1 



Aire RL h 4irR 2 



+ 



orb 



p4 
^orb 



U(R) — m e Rorb W( ), 



(102) 



that defines a (still continuum) family of circular tra- 
jectories that are stationary if we disregard (as we have 
done) the terms causing our "secondary oscillation" (2nd- 
Osc). 



D. Rad/abs balance: Puthoff's condition 



and reordering by degree of the term in R orb , 



R orb — 1" JZ h e JZ Q^(R) = 0) 







m e 
o 2 

9 -i? 2 



47re 



orb 



fi 2 Rorb 



4tt 4tt 

+ qfi e + 3m e fi e Q 2 UJ( R ) = 0, 
(107) 



where again we recall jl e — . 

A very important remark. For the system of equations 
(I104p and (|107l) (or any of the latter versions) there are 
now only two independent variables: cjir) and R or b (we 
already got rid of u)( ))- 



In the former equation, there are three free variables: 
w (o)i ^(R) an( i Rorb- We have an extra condition on L0( o ) 
that we have not used. This condition is related to that 
balance of radiated/absorbed power. For a stationary 
orbit, then, we have the condition that the emitted and 
absorbed radiating power must agree, when summed over 
the whole set of points of the orbit. A similar approach, 
as we already pointed out, was adopted previously by 
Puthoff [2l[: "It is hypothesized that (at the level of 
Bohr theory) the ground-state orbit is a ZPF determined 
state, determined by a balance between radiation emitted 
due to acceleration of the electron and radiation absorbed 
from the zero-point background" . Therefore, 



(Prad) (orbit) — {Pab) (orbit) 1 



(103) 



from where, using the expressions in ([75 ]) -([5D ]) . finally, 
Puthoff arrived to the condition: 



m e U( )R 2 orb 



(104) 



with m e the mass of the electron, u)i Q \ the orbital angu- 
lar velocity and obviously R or b the orbital radius. This 
equation obviously quantizes the (ground state) value of 
(each of the projections of) AM. 



E. An stationary trajectory is singled out rad/abs 



F. (Recovering) the Bohr radius 

We include the following calculation to add, at this 
point, clarity and coherence to the full development. 
Therefore the reader can skip it and go straight to the 
next subsection. In any case, we are not claiming here 
that the Bohr radius has anything to do with the actual 
ground state of the hydrogen. Indeed, later we will see 
there are different kinds of stable trajectories with van- 
ishing angular momentum (AM), that are clearly much 
feasible candidates for that (see Sec. \X\ for example). Nor 
is our model, in any way, equivalent to the well known 
model of Bohr. On equation (|106|) . the fist term on the 
left hand side (electrostatic force) and the lone one on the 
right (centripetal term), are clearly dominant, by several 
orders of magnitude. Equating then we have 



47re 
and therefore 



q2 R 2 = 

^orb 



h 2 R orb q 2 _ H 



2 h 2 

. Rorb=—, (108) 
47re m e 



R 



iweoti 2 

orb o 

q m e 



(109) 



which is nothing but the Bohr radius for the hydrogen 
atom. 



Puthoff's condition can be substituted in equation 
or (|102|) . yielding: 



4^o R 2 orb ^R 4 orb R orb 



+ ~54 — W (R) 



™.Rl rb 



(105) 



and, working with the common factor R^ rb , we have 



47T6, 



-±h nrb ~r 



47T 



3— Q 2 u 



h 2 R D 



(R) 



q 2 d2 , - , o - r\ h 2 Rorb 
R orb + We + 3/i e <y 2 W (i ?) = . 



4ire 



(106) 



G. Particle structure and uj 



We are concerned with the stationarity of trajectories 
with a certain stationary value (^(r) = W(^)| st . Once 
again, it is time to recall our definition of our as W(r) = 
uj — W( ), as well as note (34J. Now, for this stationary 
behavior, it is natural to assume oj\ st — 0, because the 
mass of the proton (the charge distribution) is almost 
infinite in relation to the mass of the electron. With this 
choice, 



J (o)\ 



(110) 



Some comments regarding the model we choose for the 
structure of the particle are due here. At least for the 
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simplest model of a uniform charge distribution that ro- 
tates, u = would imply that there is no net magnetic 
moment. Nevertheless, that was only a model simple 
enough to serve our purposes. From now on, we will 
adhere to the quark model: quarks are the carriers of 
the overwhelming proportion of the magnetic moment 
of the particle, the contribution of our "residual" spin 
being sufficiently marginal. In other words, the quarks 
that compose the proton add no net "inner orbital move- 
ment" to the overall spin (inner AM) of the particle. On 
our classical model, this must only happen "on average" : 
that residual classical AM of the proton may oscillate 
back and forth, but its mean value will stay zero. 



H. Higher order approx. to the simplest stationary 

orbit 



Again, we include a calculation (in this case a refine- 
ment of the former one, in Sec. HXFp that is not strictly 
necessary for our development. Anyway, we wanted to 
insist on the fact that our work stems from the inclusion 
of high order terms from the multipole development, usu- 
ally not included in the conventional models of the atom. 
Now, we do not drop our extra term in the stationar- 
ity equation. In the former section we worked out the 



condition 0J(r) 



st 



->(o)\ s ■ This simple condition im- 



mediately discretizes the spectra: up to now, we had a 
system of two equations for three variables (R or b, 
and but now those variables are only two (R or b 

and W( )). 

It is clear then that our stationarity equation (or sys- 
tem of two equations) can only be satisfied by one value 
of the pair R or b, W( a \, for the ground state, and later for 
each one of the excited states, which we will treat later. 
From (fTU7|> . with cj {R) \ st = -u 



( y- 



1 



p2 



+ qfM e - 3TO e /} e Q 2 W( ) = 0, 



I. Stability criteria for the simplest orbit 

1. Stability of dynamical systems 

We have already shown that there is a stationary tra- 
jectory for our dynamical equations, but we still do not 
know if this accounts for an "stable" one. Stability im- 
plies stationarity, but the opposite is not true. A sta- 
ble trajectory must show, aside from being stationary, 
some "resistance" to be driven apart at least against 
small perturbations. Of course this last concept of sta- 
bility only make sense when the dynamical equations of 
the system (up to now, strictly deterministic) are supple- 
mented with some stochastic terms, representing the net 
difference between radiated and absorbed power. These 
stochastic terms introduce both the possibility of energy 
dissipation or absorption. They were already included 
in Sec. IVI B 51 Of course, they only intend to represent 
roughly this processes of rad/abs from the background. 
Much more detailed calculations could be done (for ex- 
ample, we can express instantaneous radiated power as a 
function of acceleration), but they are not necessary for 
our main purpose here: we want to show that stability is 
present in this simplified (and more general) picture. 

To study stability, it is convenient here to define a re- 
duced configuration space (we are not interested in vari- 
ables such as position or angle or rotation, but only on 
the minimum set of then that characterize an orbit). On 
that (reduced) configuration space, a vector x gives the 
state of the dynamical system for a given time t. We 
thus define the following vector of dynamical variables, 
expressing the (instantaneous) state of the system: 



Rorb 

(o) 
V r 

0) 



(112) 



(111) 



where, depending on our purpose, we could omit or add 
dynamical variables. Now, given the equations of move- 
ment remain invariant with time, we can write the' 'au- 
tonomous" system of equations: 



Later we will provide a strict condition of stability of 
this first circular "attractor", that, nevertheless, is not 
yet going to correspond to the actual lowest energy state 
(a p-state, because there is net AM), mainly because we 
have ignored all the way our 2nd-FL that would introduce 
an oscillation in v r °^ . 

But we have to address some other previous questions. 
First, we have to prove that stability is (at least) mathe- 
matically possible. Later, another of this questions would 
be to open room for the explanation of the whole spec- 
trum of s-orbitals, through a different type of "station- 
ary" trajectory that must bear no average (orbital) AM. 



x = /(x), x £ R p , f : RP^R p . (113) 

Given that the function / may be (indeed it is) non linear, 
we can linearize it around a point Xo, 

x = A/,-(x-x ), (114) 

where Mi is a matrix of constant coefficients resulting 
from the liberalization. 

A necessary and sufficient condition for a stable 
"point", i.e., a circular orbit in our formulation, is sim- 
ply that the p eigenvalues of Mi must be negative. 
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2. An ideal circular orbit 



Around a certain point R or b — R^l-h we wn ^ define the 



following state vector: 



With k r , kt,k u ,A, B,C,D > 0, all roots have a negative 
real part. This is a necessary and sufficient condition for 
the point x st to be a (local) attractor in the configuration 
space of the (linearized) system. 



T)St 



,(°)\st 



(R)\ 



(115) 



with vf\ st > and, for the case already treated in 
SecIXGl u; {R) \ st = -w (o) | st . 

We can now linearize the dynamical equations of the 
system given in ([90 ]) -(|88 | and (|9Tj). Following for ex- 
ample ([53]), close to R orb = Rf rb we have {Xr(R or b)) ~ 
—k r {R or b — i?o* b ), k > 0, where the mean value stands 
for an average on an interval, but small enough to com- 
prise just a reasonable small segment of the (cyclic) orbit. 
Therefore, our equations, following now, are valid "on an 
average" along a short interval: 



J. Absence from the quantum spectrum: 2nd-Osc 

Again we stress that this lowest circular orbit is unreal- 
istic, due to the fact that we have ignored the secondary 
feedback loop (2nd-FL) that would introduce necessar- 
ily oscillations in Vr°\ Later we establish a correspon- 
dence between a set of stable orbits and the actual quan- 
tum orbitals, defining a labeling using certain integers 
n Sl n p = 0, 1, . . . that would correspond to the principal 
quantum number n. The case that we have just analyzed 
here would correspond to n p = 0. 



X. STABILITY WITH VANISHING 
(AVERAGE) AM: TOWARDS THE S-ORBITALS 
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defining the coefficients 



A = 

B = 

C = 

D = 

E = 



-qBi + 
3/t e <3 2 



2K (o) ) 



R 



p4 
n orb 



■x st ' 



3/i e Q2 



-kR orb - qEi - m ( 



K (o) ) 2 

Rorb 



(117) 

(118) 
(119) 

(120) 
(121) 



where, as indicated, each coefficient is evaluated in the 
point of stability (the point around which we are lin- 
earizing), and therefore, E\ = Ei^R^), etc. We must 
also note that all quantities are positive defined, i.e., 
k r ,kt,k u ,A,B,C,D > 0, except for E. The eigenvalue 
equation for the coefficients matrix is, therefore: 











B 



A -1 

k r A -A 

-C X + k t 

-D A + k. 



(122) 



A. Overview 

So far we have dealt with circular orbits. Clearly, these 
bear net (average) AM, in contrast to the quantum me- 
chanical lowest energy orbitals or s-orbitals. There can be 
^toTLnd, indeed, stationary (and stable) orbits (closed tra- 
jectories) with vanishing AM, but it seemed much clearer 
to us to invert the presentation (addressing first the more 
intuitive circular orbits) in the way we have done. As an 
aside, the more or less clear relation between our classi- 
cal orbits (or set of orbits) and the quantum mechanical 
orbitals will be treated in detail in Sec. IXVI Bl 

Basically, we could formulate a qualitative argument 
in the way that we did before, for the circular orbits. In 
this case the situation, however, seems a bit more com- 
plicated: this time uj [3J] does not just oscillate slightly 
in response to external perturbations, but does itself de- 
scribe an stationary oscillatory curve with a certain am- 
plitude (again, nevertheless, the value w = 0). The am- 
plitude of this oscillation should be determined working 
on the dynamical equations of the system. For now we 
will be content to say that, assuming rather natural mod- 
ifications of the "base" trajectory, we can presume that 
this amplitude may result in a reasonable value. For ex- 
ample, we can assume that the tangential velocity also 
oscillates leading to a trajectory that resembles (in the 
simplest configuration that we can consider) an "eight" , 
a trajectory whose average AM still vanishes. It is also 
important to say that to understand this kind of trajec- 
tory we need to take into account the Lorentz interaction 
on a point particle (our stabilizing mechanism can only 
produce forces in the radial direction). 

A key question is that the 2nd-Osc is again ignored, 
its role being crucial in our further justification of the 
full spectrum. To conclude, again we leave any energy 
estimation for future sections. 
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A "model" trajectory 



D. Stability: attraction dynamics 



Harmonic oscillators are not (neither classically) eigen- 
functions of an —type potential. Nevertheless, consid- 
ering the trajectories we propose are not strictly radial 
(straight lines) at all, but they combine radial and tan- 
gential components, we will adopt the harmonic formu- 
lation, with the intention of showing that (i) they are a 
"feasible" approximation to the real trajectories, (ii) they 
are stationary, (iii) a balance of rad/abs can be attained, 
hence, they can be stable. 

That last property (iii) will single out, amongst a con- 
tinuum of possible stationary trajectories, a particular 
one, the stable one or "attractor" . We will now name ujq 
the basic frequency of oscillation. This is no longer an 
angular velocity, as it was for p-type trajectories, and, 
therefore, it has nothing to do with 0J( o )- Nevertheless, 
ut \ is still perfectly defined as one of the dynamical vari- 
ables of the system, and indeed it will still be useful to 
show that stability is possible. In similarity with our pre- 
vious treatment of stability for a circular trajectory, we 
will again ignore the secondary feedback loop (2nd-FL). 
We propose a set of stationary trajectories, parameter- 
ized by an amplitude v^: 

v r° } = V r,0 COs(w i), 

and, in terms of the radial coordinate r = |ri|, 

T = Rmax COs(w i), 



(123) 



(124) 



with R ri 



.,(») 



Let us leave aside a realistic estimation of the ampli- 

( s) 

tude v^q. Instead of that, let us also suppose that, for a 
certain values of that amplitude, a power balance can be 
attained. Moreover, what we will do here is prove that 
the existence of that balance is itself enough for stability, 
given the equations that we have for the system. We will 
do this by a very general mathematical analysis. The 
main point is that, although we have v[°^ ^ 0, we do 

have (v^ )r or bit) = 0. This oscillation causes a dissipa- 
tion/absorption game but is not any more dependent on 
a particular value of v\ . Therefore, we make kt = 0. 

Besides, no "attraction" dynamics can be supposed for 
Rorb, Vr°^ or u)r R \ . Indeed, what we seek for is an os- 
cillatory one. Once the system is linearized, this kind 
of behavior corresponds to roots on the imaginary axis. 
Nevertheless, if those oscillations of R or b, v r °^ and W(r) 
are not very big (or reasonably small), we can linearize 
around the following point this time: 



Rorb 


= 


0) 
V r 


= 


0) 
V} 


= 


- 


= 



(125) 



First, stationarity is established. For that, we identify 
imaginary poles in absence of rad/abs. Therefore we do 
k r = 0, kt = and k u = 0. Our linearized system 
is similar to that of (|116[) . and the matrix of linearized 
coefficients M; such that x = M; • (x — x' st ) is written as 



Mi 



C. Stationarity 



10 

-A 1 -B' 

C 

D' 



(126) 



A key point is to see why we cannot have vf' — 0. This 
is due to the term involving u> in the equation for vf \ at 



equation (1731 or (|90|) . We will not extend much on this 
but to say that, in this kind of quasi-radial trajectory, 
the nucleus is angularly accelerated back and forth due 
to the coupling between v r °^ and oj in the corresponding 
equations. 

Besides, to understand such a trajectory, that certainly 
goes beyond the purely radial one that we should expect 
from a simple model, one has to think in the Lorentz 
interaction. This Lorentz forces go normal to the instan- 
taneous velocity, and can therefore explain (although we 
still have not calculated with what amplitude) an oscilla- 
tion such as the one we are suggesting for v\°\ Also, 
as said before, harmonic behavior is of course just a 
very rough approximation, but for the moment it will 
be enough for our purposes. 



and A', B' , C, D' obey to the same expressions as Q117p - 
(|120[) but with R or b = this time. Obviously, this as- 
sumes that the l/R or b dependence disappears for R or b 
sufficiently small, something that does occur when we 
include in the overall model the inner structure of the 
nucleus, as we are doing here. 

Again with A' , B' , C" , D' > 0, all roots A.; of the eigen- 
value problem \Mi — AI| = are purely imaginary, a 
proof of what is, once more, trivial and therefore not 
needed here. That corresponds to the oscillatory move- 
ment we were seeking for. In the same way that we had 
for our p-orbits, in absence of dissipation/absorption, a 
perturbation sends the system from one stationary orbit 
to another, both belonging to a continuum of stationary 
trajectories. Once we have established the fact that there 
is such a continuum of stationary trajectories, the ele- 
ment of dissipation and absorption from the background 
introduces the attraction dynamics. 
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XI. TOWARDS EXCITED STATES 

We have to open space for an infinite, discrete spec- 
trum, and this is done by the introduction of the sec- 
ondary oscillation (2nd-Osc), a natural consequence of 
the equations of the system, in our analysis. 

Showing that that our framework here leaves reason- 
able space to account for the existence of excited states is 
one of our main goals here. In principle, the solution re- 
duces to find new trajectories, first stationary, then nec- 
essarily stable, too, whose (mean) energy is higher than 
the already seen. So far, we have only proven the exis- 
tence of one stationary circular orbit (a candidate to gen- 
erate, in 3D, a p-orbital) and another one with vanishing 
AM (a candidate for the ground state of the system, an 
s-orbital). This has been done disregarding the oscilla- 
tion in vi°^ coming from the 2nd-FL that we have so far 
ignored. 

Precisely that 2nd-Osc is the additional element that 
will allow for an extra freedom in the power balance equa- 
tion, yielding now an infinite (discrete) set of closed tra- 
jectories that can attain that balance. Now, we intend to 
prove that excited states, with different energies, can be 
built following the same basic ideas: first we would find a 
stationary trajectory consistent with the equations of the 
system, and then we would have to prove that they are 
indeed stable against stochastic perturbations. All we do 
here is propose some particular form for those stationary 
orbits, and make some preliminary calculations on their 
energy difference with respect to the former ones. Again 
we follow, because we think is more intuitive that way, 
the same structure: first circular orbits, then the s-ones. 

The trajectories that we present in the following are 
not stable, nor are they stationary, if we stick to the sim- 
plified problem of an scalar central potential. For the 
trajectories we propose here to be possible, we must add 
both the action of the Lorentz force on a point parti- 
cle, as well as our higher order terms arising from the 
multipole orders that we wanted to include in the calcu- 
lation (our "stabilizing mechanism"). We provide here 
two sets of stable trajectories, parameterized by two in- 
teger numbers n Sl n p , ranging in principle (in principle) 
from to infinity. The notation is chosen so as to suggest 
that each group actually corresponds (all the necessary 
generalizations done: extension to 3D, stochasticity) to 
the quantum spectrum of s and p-orbitals. Indeed, as 
we will see, those integers are directly related to a cer- 
tain "matching condition" between the principal orbital 
frequency of the trajectory and the frequency of 2nd-Osc. 

Later, in Sec lXVll we will establish the correspondence 
with the quantum number n or "energy level" . Indeed, 
the correspondence is simply n s ,n p — n, though it also 
true that not all proposed values of n s and n p can give 
rise to real orbitals. For example, we must first justify the 
exclusion of the n = solutions (in principle possible in 
our scheme), as well as the one with n p = 1, as we know 
that for non-vanishing angular momentum (in quantum 
terms I > 0) we cannot have an orbital with n = 1. This 



questions will not be addressed yet in this section, though 
we will do in former ones. 

In any case, the path to prove stability of each of our 
proposed trajectories is this: First, prove stationarity 
with respect to the dynamical equations of the system, in 
absence of rad/ABS (as we have seen, for our model, this 
stationarity is not satisfied by just one trajectory but by 
a continuum set of them). As a second step, we must 
evaluate the balance between looses and absorption from 
the background: this additional condition will single out 
one trajectory from each former set of trajectories. As 
we said, the correspondence with a real quantum orbital 
is left for Sec. IXVll 



A. Parametrization of stationary orbits 

From now on wc will do, for convenience, the following 
simplification, describing both the sets of possible s and 
p-orbits with the same parameters: 

(i) From now on, R = \tq\, representing the position 
vector for the electron respect to the center of mass of the 
nucleus. The definition of our two references frames, RFq 
and RF\ , after our last relocation of mass, can be recalled 
in Sec. IV Al both RF and RFi anchored to the center 
of mass of the nucleus, but only RFq retains inertiality, 
with xi always in direction that joins that center of mass 
with the position of the electron (R or b = i? or (,Xi). 

(ii) R P ^ and lo^ will be the orbital radius and orbital 
frequency for the p-set of trajectories, therefore the am- 
plitude and frequency for the primary p-oscillation. Two 
other parameters, wjy and will describe the 2nd- 
Osc. Therefore, for p-orbits, omitting superscripts for 
clarity, 

r « [ i?o cos(wo0 + Ri cos{uj\t) } xi, (127) 
|r | « R +R 1 cos{u 1 t), (128) 

which corresponds to our expected radial oscillation. 

(iii) R^ and luq^ will be the amplitude and frequency 
of the orbital movement (primary oscillation or lst-Osc) 

(s) (s) 

for the pendulum-like or s-orbits. The pair R\°' and 
will describe the 2nd-Osc; for s-orbits, omitting super- 
scripts for clarity, 

r ~ iio cos(o;o<) xi + i?i cos(wit) yi, (129) 
|r | « [R 2 + Rl cos 2 {^t) (130) 

which corresponds to an expected secondary oscillation 
in the tangential direction, given always by yi. 

(iv) We complete the parametrization by renaming 

U(R)\max £ { w h,0! w i£o} depending on the type of tra- 
jectory. 

If no ambiguity is present, we will always drop super- 
scripts (s),(p) for clarity. We will only refer, henceforth, 
to Rq, ujq, R\, loi and uiRfi. 
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The harmonic oscillator approximation: 
decoupling of 1st and 2nd-Osc 



E. 2nd-Osc in the stationary orbits 



Talk about Coulomb well, approximate linearization, 
etc. Harmonic oscillators are classical eigenstates of the 
system, once this is linearized. Therefore, it is only an 
approximation. Indeed, that decoupling only takes place 
under that linearization, due to the imaginary poles of 
the linearized matrix. 



C. A new continuum of stationary orbits 

Therefore, in principle a stationary trajectory or "or- 
bit" is given by a quartet of parameters { Rq, uq, Ri , uji }, 
and a fifth wj^o that is (as we will see) determined by the 
former. Though we are, as announced, omitting super- 
scripts (s) and (p), it should be understood that these 
parameters do have different meaning depending if it is 
an s or p-orbit, but for our purpose here we are not in- 
terested in making any distinction. 

Now, we will impose stationarity on the equations of 
the system. In the firs place, if we ignore the secondary 
oscillation (2nd-Osc), what we obtain is a relation on or- 
bital pair (Rq,u)o): we have ojq from Rq, for instance. On 
the second place, if we do regard the secondary feedback 
loop (2nd-FL), we have is an additional relation that de- 
termines loi from i?o- This accounts for a continuum of 
orbits parameterized by Rq. However, R\ is still a free 
parameter. We bear in mind that. 



D. Some necessary context: rad/abs in absence of 
2nd-Osc 



We will not repeat PuthofFs calculations here, but just 
remark some questions of interest. For circular orbits 
with a higher R 0l the dissipation decreases, as we can 
see from this calculation... Indeed, applying (|B3[) to the 
radiated power, 



We know, because of the coupling between equations 
and (pTTj) (the coupling is also present when we do 
not include the radiative corrections), any circular sta- 
tionary trajectory is not compatible with the condition 
v r °^ = 0. Nevertheless, it is interesting to do some pre- 
liminary study when we do impose v r °^ = 0, which itself 
will help understand how excited (higher energy) states 
can be identified with the corresponding stationary tra- 
jectories. Besides, for a pendulum-like trajectory, the 
coupling between (1551) and ([TO]) forbids a stationary tra- 
jectory with v[°^ — 0. Both conditions will add a certain 

natural frequency uj\ and lo^ to the orbital movement. 
As we have said, it is this 2nd-Osc that will make possi- 
ble to deviate from PuthofFs unique balance condition, 
allowing for an infinite set of stable orbits, discretized by 
a certain matching condition. 

Again, later, we will establish a correspondence be- 
tween a labeling number n (similar to the principal quan- 
tum number) and each excited trajectory, and it is no 
less important to remark that the presence of this 2nd- 
Osc will explain why no p or s orbits exists for n = 0. 
Moreover, basing our argument in geometric considera- 
tions, we will exclude the p-orbit for n = 1, too, therefore 
in coherence with the actual atomic spectrum. 

PuthofFs condition as given in j21| only allows for one 
unique stable circular trajectory (and one unique pen- 
dulum one). The only possibility to enhance the set of 
stable trajectories demands the introduction of new de- 
grees of freedom in the orbit: that way, some degrees of 
freedom can radiate more or less than they absorb, coun- 
teracting the excessive or defective absorption of the oth- 
ers. This said, we have to realize that it is the structure of 
the system itself (the dynamical equations) that provides 
us with those extra degrees of freedom. As we have seen, 
there is a coupling between equations (|88l) and (|9Tj). so 

trajectories of the kind Vr — can never be stationary 
at all if 7^ 0. 



Rq 



2 



{Prad) 



1 



47T£o 
„2 p2, ,4 

67reoc 3 



2 ■ 





rrieOJ, 
1 



m e Rl ' 



r?4 : 



(131) 
(132) 
(133) 



which means that for feasible excited states correspond- 
ing to circular orbits, the dissipation must be compen- 
sated by a significant contribution from the secondary 
oscillation, that in turn will not cause much absorption. 
A similar conclusion can be obtained for pendulum-type 
orbits. 



Besides, the coupling between (155)) and (|50|) . that 
stays irrelevant in stationary circular trajectories because 
= 0, acquires key importance in trajectories of the 
pendulum type, making impossible that absence of ra- 
dial velocity any more. To the orbital movement (either 
if it is circular or pendulum-like) we have to add an os- 
cillation on Vr . An obvious candidate is a harmonic 
oscillation. Why? Because the (linearized) system has 
imaginary poles in its frequency representation (extend 
on this), and harmonic oscillation therefore is indeed an 
"eigenstate" (we mean a classical one, here) of the sys- 
tem. In the following, let oj\ be either Wosc or wise, i.e., 
the 2nd-Osc frequency respectively for the s and p-types 
of stationary orbits. 
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F. Necessity of a matching condition 

In principle, stationarity requires a "matching" condi- 
tion: oj = u>i/N for some integer N = 1, 2 Later, we 

will try to establish a correspondence between the integer 
N and the well known principal quantum number. 

G. Feasible p-orbits 

We provide now an expression for the whole dis- 
crete p-spectra, introducing now, as already advanced, 
a parametrization by a certain integer n p , that later we 
will make correspond to the well known principal quan- 
tum number. Initially, 

Vt a) ~ U( )R orb , (134) 

» v^cos^t). (135) 

1. p-matching condition 

Now, our matching condition will be 

ui = n p uj {o) , n p = 2, 3, . . . , (136) 

where the absence of n p — 1 corresponds to the exclusion 
of the p-orbitals for the first atomic level (later we will 
give some more justification on this) , yielding 

ps uj( o) R orb , (137) 

vi o) w v^cos(n p u (o) t), (138) 

where, clearly, the higher the n p , clearly the higher the 
energy. Besides, the amplitude of secondary oscillation 
will possibly depend on n p , determined by the dynamical 
equations of the system: 

«3 = 2,3,... (139) 




FIG. 5: In our proposed p-orbits, the orbital distance 
undergoes a "secondary" oscillation (2nd-Osc) in the way 
R = Ro + Ri cos(cjit) (which is simultaneous with an os- 
cillation in the AM (or angular velocity) cj of the nucleus, 
vanishing on average). It may look a counterintuitive 
trajectory for a central potential, but mind the in- 
tervention of the "inner" degree of freedom! The fre- 
quency matching condition for a p-orbital is u\ = n v u)o. This 
basic p-orbit gives rise to a rotationally invariant electronic 
density around the axis zi. When we analyze the situation in 
3D, we have to consider a modulation of our stabilizing inter- 
action from the misalignment of that axis zi and the magnetic 
moment of the electron [i e (this misalignment is nothing but 
the classical counterpart of the quantum LS term). 

H. Feasible s-orbits 

In the way that we have done for the p-spectrum, we 
introduce a parameter n s , which will later allow us to 
establish a correspondence with the well known principal 
quantum number. Initially, 

R i? cos(w t), (140) 
v[ o) « ^gcosM). (141) 

I. s-matching condition 



2. Stationarity and stability 

To proof stationarity, it is enough to show that the 
system, in absence of rad/abs terms that will further on 
privilege a discrete set of orbits and explain their attrac- 
tive behavior, has a frequency description where all poles 
are purely imaginary. To analyze it, we first had to lin- 
earize around a particular point of equilibrium, but once 
the secondary oscillation dynamics is included, neither of 
the dynamical variables used in our description of the sys- 
tem (see...) adopt a constant value around the stationary 
orbit. Anyway, we can keep on considering that the oscil- 
lations in v[°^ are small, and the description already used 
in Sec lIXI still remains adequate. The attraction behavior 
is again provided once we include the rad/abs terms. We 
would need to prove that this new balance is attainable. 
We do not consider this a major difficulty, however, and 
we leave any calculations for the Appendix. 



Now we have 

wi = 2n s uj Q , n s = l,2, (142) 

yielding 

R « i? cos(w t), (143) 

^ (o) pa v t ( * ) coB(2n.w ) ) (144) 

where, clearly, the higher the n s the higher the energy, 
and where the amplitude of the 2nd-Osc, this time affect- 
ing what we defined as the tangential velocity component, 
will possibly depend on n s , determined by the dynamical 
equations of the system: 

4o=^oK), n s = l,2,..., (145) 

Special attention should be taken to the factor 2: we 
need a trajectory that "crosses " with itself, so the aver- 
age AM over a whole cycle yields exactly zero. Again, 
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the reason for the absence of solution n s = is that it 
disregards the 2nd-Osc so typical of our system. Follow- 
ing our program, as we did for the p-case, we now have to 
check if they are indeed stationary, and then if they can 
be stable (strict stability happens whenever a power bal- 
ance can be attained, singling out a particular trajectory 
from the set of all stationary ones). 




FIG. 6: In the proposed s-orbits, the secondary oscillation 
(2nd-Osc) takes place in yi, as well in small variations of 
the angular momentum (AM) of the nucleus (that vanish, in 
average). The average overall AM of this orbit is zero as 
expected for the quantum mechanical counterpart. The fre- 
quency matching condition for an s-orbital is uii = 2n 3 uuo, 
with n s — 4 in this example. It may look a counterin- 
tuitive trajectory for a central potential, but mind 
the intervention of the "inner" degree of freedom! 
Of course this trajectory is not stationary in the purely elec- 
trostatic picture, but it is feasible when we include both the 
Lorentz and higher order terms. On the other hand, this basic 
s-orbit gives rise to a rotationally invariant electronic density 
in 3D, when the axis xi is allowed to precess freely in response 
to stochastic perturbations. 



2. Stationarity and stability 

Our reasoning is the same as for the p-spectrum: if 
we consider that the secondary oscillation (2nd-Osc) is 
small enough, then the linearized description of the sys- 
tem used in Sec. IXl remains approximately correct, and 
the purely imaginary character of the poles in the fre- 
quency representation, hence the purely oscillatory be- 
havior, guaranteed. Once the existence of a continuum 
of stationary orbits is proved for the equations of the 
system in absence of dissipation/absorption terms, their 
inclusion privileging a discrete set of orbits (and explain- 
ing the attraction behavior). 

Our next task will be to do a first approach to the 
energy spectrum. We have already given the basic equa- 
tions for the stable orbits (for the linearized system, in 
2D and etc.), parameterized them by a natural number 
n s or n p , and now we want to estimate the quantity 
AE(n a ,n a + 1) = E(n a + 1) — E(n a ). I.e., we want 
to estimate the energy gap between adjacent states. 



XII. A FIRST APPROACH TO THE 
E-SPECTRUM 

In principle, we would have to calculate the following 
energies, for p and s-orbits, respectively, 



[P \n p ) + E^{ 
E^(n s )=E[ s \n s )+E ( 2 s \n s ), 



(146) 
(147) 



where the subscripts 1, 2 make reference to the principal 
(orbital) or secondary oscillatory movements that char- 
acterize each orbit. We will hereon omit superscripts for 
simplicity (the dependence on n s or n p eliminates any 
ambiguity) . 



A. The p-spectrum (in absence of 2nd-Osc) 

For p-orbits, disregarding the 2nd-Osc, we know from 

E3), 



Ei{n p ) = -E k in(n p ) = -^m e (v^) 2 = -im e cj 2 i?g, 

(148) 

and now, applying the frequency matching condition, 

Ull = TlpLVo, 



1 ui 2 
Ei{n p ) = --m e -^Rl, 



but, from ([B28| . 



so, therefore 



Cl>i OC 



1 1 



« „2 R 2 



Considering that 



1 1 1 

"52 ~ "52 - 2 "53 ( R -Ro) + -- 
n tl n 



(149) 



(150) 



(151) 



(152) 



with R(n p ) - R(n p = 2) << R(n p = 2) 3 , Vn p > 2, we 
can make the approximation 



1 



1 



Ro(n P ) 2 R {n p = 2y 
and finally arrive to 



n p > 2, (153) 



ASi(n p ,n p + 1) = Ex(n p + 1) - Ex(n p ) 
1 1 



n% (rip + 1)2 
and therefore AEi(n p , n p + 1) > as expected 



, (154) 
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Doing a finer estimate, however, from Appendix \B_ 
we can see that actually 



1 f 



Ei(n p ) oc -^y- 



(155) 



but one has to bear in mind that all these estimates are 
done on a two-dimensional model of the problem, so their 
significance is only relative. 



B. The s-spectrum (in absence of 2nd-Osc) 

For s-orbits, disregarding the 2nd-Osc, we know from 
(1B161) . 

E x {n s ) = -E kin {R = 0) = -~m e oj%R%, (156) 

but applying the frequency matching condition, u>i — 

n s oj Q , 



£i(n P ) = \ m ^ R l + E vot{R = 0), 



but, from (fB28|) . 
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so, therefore 
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Considering that 
1 1 



Rq Rq i?n 



(R - Ro) 



(157) 



(158) 



(159) 



(160) 



with R(n s ) - R(n s = 1) << R(n s = l) 2 , Vn s > 1, we 
can make the approximation 



1 



R (n s y R (n s = iy 
and finally arrive to 



n p > 1, 



(161) 



A£a(n s ,n s + 1) = Ex{n s + l)-Ex{n 
1 1 



(162) 



Actually, if we do a finer estimate using the results of 
Avvendix \B 31 we can see that actually 



1 

Ro 



OC Us 



Ei(n s ) oc |-, 



(163) 



where again we remind the reader that all these estimates 
are done on a two-dimensional model of the problem, and 
under many simplifications, so their significance is only 
relative. 



C. 2nd-Osc in the E-spectrum? 



For the moment, we will not include E^ (%>)> ^2"' ( n s) 
in the calculation. There is a reason for this. Our 
secondary oscillation (2nd-Osc) expresses a resonance 
behavior: the system has the capacity of absorb- 
ing/dissipating whatever power is necessary to keep its 
orbital movement at a particular frequency (wi). There- 
fore, the energy corresponding to that 2nd-Osc is, always, 
"immediately" absorbed or given to the background, and 
with this we mean that in a much smaller time scale than 
the one corresponding to the principal oscillation (that, 
from the point of view of this 2nd-Osc, exhibits a practi- 
cally constant behavior) . This is the main reason why we 
believe this energy should not be included in the estimate 
of the energy spectrum. 



XIII. NON CIRCULAR ORBITS IN 2D: 
ELLIPTICAL. ORBITS WITH HIGHER AM 



We have so far only addressed the (trivial) case of cir- 
cular orbits in the plane (2-dimensional problem). Bohr 
used elliptical orbits to account for states with quan- 
tum number I > 1. This would be the starting point 
in our approach, as well. More complicated stable orbits 
( "limit cycles" in the literature of systems theory) , given 
by the nonlinear nature of the dynamical system, require 
a more sophisticated analysis, but some simple "mind ex- 
periment" (giving some initial condition vi°^ ^ in the 
planar problem, for instance) may surely convince the 
reader that these do actually exist. 

Of course, once we have a nonlinear behavior, stability 
is restricted to certain ranges for initial conditions. This 
is another task. Comment on the possibility of elliptical 
attractors. Use of the Poincare method. 



XIV. A 3D MODEL 

Only some general, preliminary ideas here. Convenient 
to include it anyway. 



A. A first approach to p-orbits in 3D 

When we consider general initial conditions, the new 
forces and/or torques (to the leading order) are strictly 
normal to the ones in the planar case, so the stability 
mechanism remains valid. To make our picture consis- 
tent, the classical counterpart of the quantum LS inter- 
action must appear here as one of the new terms. This 
interaction would modulate the strength of our stabiliz- 
ing mechanism: it decreases its intensity if the magnetic 
moment fi e is not aligned with the axis Zi for the circu- 
lar orbit. On the other hand, we know that a dipoling 
interaction between magnetic moments of electron and 
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proton must be also present (although it does not ap- 
pear as a primary term in the Hamiltonian expansion), 
and it tends to anchor the relative direction of those two 
moments. The combination of this last dipoling inter- 
action with the classical counterpart of the quantum LS 
must be the key to determine the probability distribution 
of some "preliminary" p-orbitals in 3D. We say prelimi- 
nary because we are also considering all the way here that 
the inner structure of the nucleus is perfectly isotropic: 
nothing further from truth. Perhaps those differences 
may well explain the lack of rotational symmetry (not 
even around the axis of the orbit) of the "real", quan- 
tum, p-orbitals. 



B. A first approach to s-orbits in 3D 

For this kind of trajectories we directly include here 
some ideas on the extension to 3D. On physical terms, the 
situation we have is one where the electron oscillates back 
and forward, moving in a direction that, due to stochastic 
interaction with the background, can freely precess giving 
rise to an spherically symmetric set of trajectories. Of 
course, reorientation of the spin of the electron is needed, 
but this is the consequence of the terms in Sec. II V Jl This 
is the kind of behavior we would expect for a quantum 
s-orbital. Hay que elaborar mas la interpretation de cada 
termino en Sec. \IV J\ 



Consider an origin of coordinates Oi and let ri be the 
position vector defined in respect to it, as well as a second 
origin of coordinates 2 such that the position vector is 
now r 2 . Obviously, the following difference 

ri, 2 =r 2 -ri = Oi-0 2 , (164) 

is clearly a constant vector. Now, let A be a closed tra- 
jectory, with a period T, so r,(T) = r,(0), for i = 1,2. 
For the movement of a system of mass m whose posi- 
tion is (simultaneously) given by i"i(t) and r 2 (t), in the 
respective system of reference, it is easy to prove: 

(L 2 )a = i J L 2 dt = - J r 2 A p dt 
1 f T 

= Tj ( r i + r i,2 ) Ap dt 

1 f T 1 f T dri , 

= Tj riApdt+^A^ m— dt 

= (Li) A + yri, a Am[ri(T)-ri(0) ] 

= (Ma, (165) 

Thus, the classical AM averaged over a closed orbit is 
independent of the origin and can be added up, just like 
its quantum counterpart, which suggests a possibility of 
relating both concepts. For instance, for a p-orbit A p , 
(Li) Ap = (L 2 ) Ap = h/2, and for an s-orbit A s , (Li)a. = 
(L 2 )a s = 0. 



XV. SIMULATIONS 

Numerical simulations. No simulations for the mo- 
ment. 



XVI. THE QUANTUM HYDROGEN: ATOMIC 
ORBITALS 

We now relate our results to strictly quantum features 
of the real hydrogen atom. With this section we complete 
the bridge from our SED context to QM. 



A. Classical and quantum AM 

Classical and quantum angular momentum (AM) are 
often treated as different concepts with only a limited 
relation between them. If we let aside rotations that are 
not (globally) isomorphic to the group SO (3) of rotations 
in ordinary space, the main difference between the con- 
cepts of classical and quantum AM is that the former can 
not in general be added, because of their dependence on 
the chosen point in respect to which they are evaluated. 
However, the average value of a classical AM on a closed, 
periodic trajectory is independent of the origin, as we will 
see from the following calculation. 



B. On the concept of classical "orbital" 

Most of the time, we have been talking about trajec- 
tories, rather than orbitals. We feel in need to establish 
a clearer and more convincing bridge from one concept 
to the other. We have worked with orbits, or closed tra- 
jectories, in the plane. As a matter of fact and keeping 
an eye on the quantum orbitals, the extension to the 
3D problem arises the need to generalize that concept 
to a new one, that of an "orbital". We probably said 
before that with the term "orbital" we mean a collec- 
tion of points in the configuration space of the system, 
meeting the condition that, if the system is inside, it 
would probably stay inside with overwhelming probabil- 
ity. With this concept we are going two steps forward. 
On one side, where needed we generalize from one par- 
ticular trajectory to a set of them (for example, a set of 
s-orbits containing all possible directions, so the result- 
ing orbital exhibits spherical symmetry). On the other, 
we introduce the stochastic character; those stationary 
orbits, when subjected to the action of a perturbating 
background, gives rise to a probability distribution. 

Besides, so far we have only been always talking about 
"orbits" as "closed" or "periodic" trajectories. In regard 
to this, we must say that "stability" can also be seen as 
a consequence of a more "chaotic" dynamics, for exam- 
ple through trajectories that adopt the so-called "fractal" 
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kind of behavior: it never goes through the same point 
(initial condition or state) twice, but the evolution stays 
forever into some particular region of the configuration 
space (space of values for the dynamical or state vari- 
ables) of the system. At this point we have already shown 
that, at least in an idealized model, this purely classical 
stable orbitals exist, and they have a discrete character, 
in the sense that each of them correspond to a particular 
stationary trajectory (or a set of them). This is indeed 
the case of our s-orbitals. By contrary, p-orbitals in 2D 
arise from one particular circular trajectory. Their ex- 
tension to 3D requires to consider a modulation through 
the classical counterpart of the quantum LS effect. 

C. Equivalence to quantum orbitals 

1. The "node" problem 

(to be updated: the LS argument does not work in the 
I = case, but therefore we can invoke Fritsche's 
reasoning: pure states cannot be prepared) 

Ultimately, our classical orbitals are characterized by 
a probability distribution of the electron around the nu- 
cleus. The density of probability may extend to infinite, 
due to the fact that there is a non-vanishing probability 
that high values of the background field introduce such 
strong oscillations in the particle movement (anyway, the 
bulk of the distribution must be confined to a finite region 
of space). Therefore, our picture can find quite a consis- 
tency with the purely quantum one. Besides, it must be 
said here that, while non relativistic orbitals in the hydro- 
gen problem show disconnected regions and nodes where 
the probability density vanishes, this is not the case of the 
fully relativistic ones (for instance see ]3(\]). 

2. Phase averaging and symmetry 

(to be completed) 

For both types of stationary orbits, the whole spec- 
trum of stable trajectories we just proposed exhibits sev- 
eral free parameters: the phase (for a given principal 
axis) in the case of the circular ones, the principal direc- 
tion in the second case (s-orbits). Now, in both cases, 
when averaged over all possible angles and phases, they 
recover the symmetry of the lowest one. Their associated 
distribution of probability is centered, nevertheless, at a 
different orbital radius, as the point of balance of radi- 
ated and absorbed power is this time different. This is 
an important step towards the quantum orbitals. Now 
we recall Figs. |5] and [5j With those in mind we present 
the following two figures: [7] and [3 

The extension of the s-orbits to 3D had already been 
discussed before, in Sec. IXIVB1 About the p-orbitals in 
3D, we can only say, for the moment, that clearly we have 
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FIG. 7: The basic s-orbit that gives rise to a rotationally 
invariant electronic density in 2D, when the axis xi is allowed 
to precess freely in response to stochastic perturbations. It 
may look a counterintuitive trajectory for a central 
potential, but mind the intervention of the "inner" 
degree of freedom! 




FIG. 8: The basic p-orbit that gives rise to a rotationally 
invariant electronic density around the axis zi. When we 
analyze the situation in 3D, we have to consider a modulation 
of our stabilizing interaction from the misalignment of that 
axis zi and the magnetic moment of the electron fi e (this 
misalignment is nothing but the classical counterpart of the 
quantum LS term). 



a trajectory that is associated with a particular direction 
is space. Whether there may be only three orthogonal 
directions giving rise to simultaneous stable trajectories 
is a question that we will left aside for the moment. 
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3. n s , n p against the quantum number n 

With n the principal quantum number, the correspon- 
dence is completely straightforward, 

n p — > n, (166) 
n s — > n, (167) 

with n p and n s as defined in Sees. IXI Gl and IXI HI It is 
remarkable that in our framework, no n s — or n p = 
stationary (hence not stable either) orbits exists, due to 
the presence in the equations of the 2nd-FL that we have 
already commented in Sec. I VII Bl for instance, and there- 
fore we do not have to exclude that possibility as an "ad 
hoc" hypothesis. Another question is why should we ex- 
clude the n p — I case (the first quantum p-orbital corre- 
sponds to n = 2). A closer inspection of the equation for 
the orbit leads to think that, due to obvious the lack of 
inversion symmetry when n p = 1 , the implicit decoupling 
of the movement in an orbital (primary) and secondary 
oscillations (2nd-Osc) cannot be assumed any more. 



D. Orbital double occupancy: ideas 

A flip of the electron spin corresponds to a sign inver- 
sion for the classical magnetic moment of the electron, 
from fi e z to — /i e z. As a consequence, terms correspond- 
ing to the Lorentz force and our Inverse Magnitude Spin 
Orbit (IMSO) coupling reverse sign. This accounts for 
a shift in the energies (the electrostatic terms obviously 
remains constant), but in a way that stability remains 
equally valid, both for the s and p-case. In our conclu- 
sions we will comment more on this issue, even suggesting 
an interpretation in terms of a classical Pauli principle, 
that now arises as a consequence of our model rather than 
being an "ad-hoc" hypothesis. 



XVII. DISCUSSION 

Here, we try to add clarity on some points. We may 
also comment on alternative approaches to some of the 
questions treated so far. 



A. Rad/abs in the dynamical equations 

Perhaps we would need to dedicate some more atten- 
tion to the way in which we include rad/abs in the dy- 
namical equations of the system. This is done in equa- 
tions {SO}, (J59J and (HU, by means of the new terms dHFJ, 
([53"|) and ([57)1 . Those terms represent an unbalance of ra- 
diated and absorbed power in a quasi-instantaneous, but 
not instantaneous, basis. Therefore, they are averages, 
but over an interval small enough so that the dynamic 
included in (JOT), (EH) and (HJ) stays, effectively, "frozen". 



Looking at (|78|) - (l80l) . we see the first of them is a cubic 
polynomial on W( Q ) , while the second is a fourth order one 
(both positively defined for uj^ positive). Clearly, this 
two curves do cross each other, the point where they do so 
being modified by the extra factor R 2 orh - This obviously 
gives rise to some "attraction" around a certain value 
wfly Above this value, the probability of loss prevails. 
Below that value, the probability of absorption prevails. 
We must recall here that, in our formulation, is not 
one of the state variables of the model, but the tangential 
velocity component v[°^ — uj^R or f,, and also uj — uJ(r) + 
W( ) are so themselves. 

We are interested, at least primarily, on circular or- 
bits. Stationarity of these imposes an extra relation on 
the pair R or b, ^(o): aside from angular momentum (AM) 
quantization, already obtained from the balance of aver- 
age power around a closed orbit. This condition results 
in the corresponding (discrete) set of pairs of values, and 
therefore, it is justified to employ a prevailing dependence 
on R or t for \r, as it is done in (j53"|). This dependence can 
be linearized around the point of interest. Strictly speak- 
ing, we know that Puthoff's equations (|78|) - (l80l) only ap- 
ply to those circular orbits, but we have extrapolated 
that behavior to a situation where the system is varying 
dynamically the parameters of its orbit, and therefore 
this does not have to remain circular or stable at all (its 
stability or instability being precisely the ultimate object 
of our study here). The extrapolation seems natural in 
the case of (|53"|) . anyway. 

Less evident is our choice for ([5o| and (1571) . For the 
first of them, we can simply state that because we have, 
as defined, — W( )i? or t,, an attractor in u)fK induces 

one in v^\ st — ujfKR^. b . In the case of (l87l) . we have 
that uj — CLi( fl ) + aj( ) is none other than the angular ve- 
locity of rotation of a charge distribution around its axis, 
always in respect to an inertial frame of reference. There- 
fore, the rad/abs game applies again, and it is natural to 
assume that an attractor must exist at a certain value 
uj = uj\ st . Indeed, we have seen, for the simplest of sta- 
tionary states, uj\ st = W( D )| S *. 

The reader will agree that this is still hardly a rigorous 
treatment, but we are content to show that stability is 
not only possible, but arises as a natural behavior purely 
from the dynamical equations of the system. Though we 
have made some assumptions, they are rather natural. 



B. Relativistic considerations, "rigidity" and 
self-reaction forces 

Within our context here, all the discussion about the 
implications of a point particle in electromagnetism is 
unnecessary. On one hand, and sufficiently away, the 
electron is seen as a point-like magnetic dipole (which is 
conventional, again, sufficiently far). On the other, far 
but not infinitely far, the proton is seen as a distribution 
of charge (and magnetic momentum) that can "rotate" 
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as a whole. Besides, any reference to the rigidity of the 
distribution can be regarded as a device to make cal- 
culations simpler, but, qualitatively speaking, makes no 
difference in relation to our main results. As matter of 
fact, throughout our development in Sec. IIV1 in no mo- 
ment we have considered the forces that some elements 
of charge in the distribution exert over other elements 
in that same distribution; indeed, the mechanism we ex- 
pose is dominant up to a certain order, whenever the two 
objects remain sufficiently apart. Higher orders may de- 
form the shape of the distribution, or have other effects 
that, whenever the range of distance is the appropriate, 
do not affect our main argument. 

We have nevertheless already briefly commented on the 
(rclativistic) difficulties that arise, not only (i) when we 
deal with "rigidity" , but also (ii) with the picture of spin 
as a spinning sphere, rotating around its axis, if we try to 
adjust the angular velocity so as to reproduce the (phe- 
nomenological) magnetic moment of the proton. The dif- 
ference there is made in what kind of distribution we as- 
sume. Simple models such as those of a rotating solid 
sphere or a spherical shell of charge are ruled out; on the 
contrary the quark model, where some discrete entities 
carry charge, and more important to our purposes, all the 
magnetic momentum, seems to solve those difficulties. 



XVIII. CONCLUSIONS 

A. General conclusions 

In this paper, we have tried to reasonably establish 
a framework to support the idea that SED (i.e., classi- 
cal Maxwellian electromagnetism plus a background of 
radiation), in combination with a model of the nucleus 
where this presents some very elementary inner structure 
(specifically, a model of the proton as a spatial distribu- 
tion of charges), is perfectly capable of explaining atomic 
stability. Perhaps we should remark, before any other 
feature of our model, the following point: an atomic tran- 
sition implies an emission/absorption of a wave-packet 
with an energy h(uii — ujj), where u>i,u>j are simply the 
orbital angular frequencies of the basic stable orbits cor- 
responding to the orbitals involved. This comes as a con- 
sequence of our very basic hypothesis, rather than being 
introduced as a principle as it does in QM. 

Along with the ideas, we have done many specific cal- 
culations that, for the moment, reinforce our previous 
qualitative approach. They must, nevertheless, be re- 
garded only as preliminary and we hope they will help, 
on future developments, to check whether we can go for- 
ward and converge to more accurate estimates of quanti- 
ties such as the hydrogen atomic spectra. The ultimate 
goal of this paper is, therefore, to settle a plausible frame- 
work, not necessarily completely quantitative. 



B. The IMSO coupling 

The first main result presented in this paper is the 
identification of a classical mechanism that appears to 
be capable of stabilizing the electronic orbitals against 
(small) perturbations. This mechanism becomes exposed 
once second order moments are considered for a distribu- 
tion of charges, under the action of a point-like magnetic 
source. This interaction couples the magnitude of the an- 
gular momenta of the orbital movement, with the magni- 
tude of the inner AM carried by the charge distribution 
rotating around its axis of symmetry. 

In Sec. Ill Dl we did a quick overview of the main 
terms in the classical (also quantum mechanical) Hamil- 
tonian for the hydrogen atom. These range, as a re- 
sult of a multipole expansion for the proton, from (i) the 
purely Coulomb term, (ii) the Lorentz interaction, (iii) 
the dipoling (or in QM, spin-spin) interaction, to (iv) 
the spin-orbit interaction (or LS in QM). That said, if wc 
now include the second order (spatial extension) for the 
nucleus, we get the terms corresponding to what we have 
called IMSO coupling. The identification of this higher 
order interaction is done in two steps. Having worked 
with a previous electromechanical model, with, in prin- 
ciple, no resemblance with the actual hydrogen atom, 
we later depart from it and reinterpret our equations in 
Sec. El Its presence in a realistic system (a hydrogen 
atom) is now established. 

Another question is why we name it by "spin-orbit" . 
This is done because, certainly, it bears quite a relation to 
the conventional spin-orbit coupling (LS term in QM) : it 
is, in order, the following correction to orbital movement 
of the electron around the nucleus. Nevertheless, we must 
stress several facts here: 

(a) First, we are talking about a coupling between the 
AM of the electronic orbit and a classical inner AM of 
the proton (nucleus), not that of the electron. Through 
more or less complicated arguments, we link this classical 
AM to a residual oscillation of the proton as a composed 
object (the quarks inside being point-like entities carry- 
ing charge and spin - inner AM) , in the form of a rotation 
around its axis of symmetry. The integral of this oscilla- 
tion would vanish, so in average the value of the classical 
(inner) AM would correspond to the quantum mechan- 
ical value: from here, our choice of the terms "residual 
spin" (RS) and "Inverse Magnitude Spin Orbit coupling" 
(IMSO). 

(b) We have to insist in that hey are terms of a different 
nature. In QM, the LS term couples an effective magnetic 
field, resulting from the orbital movement of the proton 
from a reference frame attached to the electron, to the 
orientation of a magnetic dipole (that of the electron), 
not changing its magnitude. Our IMSO interaction cou- 
ples the magnitude of the AM of the orbit (given by the 
pair of variables u)r \ and R or b) to the magnitude (mag- 
nitude of the oscillations) of that "residual spin" of the 
proton ( coming from a non-zero value of u> ), not to its 
orientation, in contrast to what the well-known LS term 
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in quantum atomic physics does (there, the coupled spin 
is that of the electron). 

(c) Probably, the contribution on an average over a 
closed orbit of this new term is always vanishing (again 
in contrast to the quantum - or classical - LS shift). A 
convincing general proof of this is not within our reach 
yet. On a quasi-instantaneous scale (transparent to QM) 
this kind of effect would indeed be observable, but at 
least not by conventional atomic spectrography. 

Aside, because our high order interaction takes place 
exclusively via the Lorentz law, it is reasonable to ex- 
pect that such correction in the energy spectrum vanishes 
identically, although this is a matter still to be proven 
more rigorously. Indeed, we tried to shed some more 
light on this question with some calculations in Sec. IA II 
(Appendix). 



C. Why is it stabilizing? 

Stability of the atomic orbitals is, at least at present, 
an exclusively quantum mechanical prediction. On the 
other hand, within the SED framework and up to now, 
the ZPF background provides a mean to counteract the 
loss of energy that a charge undergoes when accelerated, 
through radiation. However, given the stochastic nature 
of the ZPF, the cancellation is just on average. Nothing, 
however, had been said yet on how could this balance be 
attained on a quasi instantaneous basis, or how the orbit 
could accommodate the instantaneous excess or defect of 
energy without modifying the pair R or b, 0J( o ) (therefore 
keeping the stationarity of the trajectory, losing it other- 
wise). In other words, it does not explain why the system 
should remain confined to a reduced portion of its pos- 
sible configuration space, that portion corresponding to 
what we know by an orbital. 

It does not explain, neither, why there should exist a 
discrete set of energies for the orbitals with a given AM. 
However, if we provide the system with the capacity of 
storing energy on the fluctuations in that classical (resid- 
ual) "spin" of one of the particles, the former picture can 
be now made to work. This classical residual spin of 
the particle must nevertheless be able to participate in an 
interplay of energy exchange with the orbital degrees of 
freedom, and in this paper we have shown that the forces 
and torques for this are present, once that one simply al- 
lows for a slightly more complex model of the situation, 
where the nucleus is permitted to have a non vanishing 
second order moment. 

Now, the instantaneous excess or defect of energy is 
stored in the proton "residual spin" , so to say, the (resid- 
ual) AM due to orbital movement of the quarks within 
the nucleon structure. Adopting the quark model is nec- 
essary because the quarks are carriers of net AM and 
magnetic moment, avoiding that way difficulties with 
special relativity (those difficulties arise when the more 
simple model of an spherical distribution of charge is as- 
sumed for the electron, and also for the proton). 



D. Why a discrete spectra? 

PuthofFs work accounted for the necessity of the quan- 
tization of AM, but now, the spectrum of parameters for 
the stable orbits is also discrete, when we impose a con- 
dition of "stationarity" over the dynamical equations of 
the system (PuthofFs argument did not take them into 
account). For instance, for a stable circular orbit, the 
possible values of both orbital radius R or b and the angu- 
lar velocity 0J( o ) are also discrete, aside from the overall 
AM. This point is of key importance, as the previous 
picture only accounted for the discretization of the spec- 
trum of AM. But yet, we could go the opposite way. If, 
in the first place, we impose stationarity what we have is 
a continuum set of stationary orbits, all of them plausible 
candidates to be stable, then, as a second step we can de- 
mand Puthoff 's power balance, and end up with only one 
orbit that would have, therefore, perfectly defined values 
for its main parameters (R or b and the angular velocity 
UJ( ))> as we M as f or H s overall average AM. 

Of course, this does not mean that real orbitals are 
reducible to neat classical orbits in that way: under the 
bath of radiation, they oscillate around that privileged 
trajectory or set of trajectories (this is the natural situa- 
tion when we formulate the three dimensional problem) . 

E. And infinite? 

We had explained the discreteness of the spectrum, but 
this accounts for only one (feasible) stable trajectory for 
each particular value of AM. The matter of how can we 
account for an infinite set of stable orbits for each of those 
values of AM was the following question to be addressed. 

So far, we had only addressed the orbital movement 
of the electron. This movement, that we had called "pri- 
mary oscillations" , could be stabilized by a feedback loop 
(primary) relating radial velocity and relative angular 
velocity in respect to the nucleus. But moreover, the 
dynamical equations of the system also include another 
feedback loop, that we had called secondary, relating 
the radial and tangential components of velocity. This 
secondary loop introduces an additional oscillation that 
takes place entwined with the principal, orbital oscilla- 
tion. Puthoff 's initial calculation only took into account 
the primary oscillation, therefore yielding a unique point 
of equilibrium of rad/abs processes. When the presence 
of the secondary oscillation (2nd-Osc) is taken into ac- 
count, an infinite set of points of equilibrium becomes 
possible: the excess or defect of power in the primary 
oscillation can be exactly canceled by the absorption or 
dissipation through the secondary movement. 

In principle, that infinite set would be continuum, as 
well, but the inclusion of an additional, and rather ob- 
vious from the point of view of the stationarity of the 
orbit, restriction relating the frequencies of both oscilla- 
tions. This restriction recovers the discreteness of the 
spectrum, but now that spectrum includes an infinite 
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number of stable (classical) states, with increasing en- 
ergies. 



F. Stable orbits with and without net AM 

The second of the main results of this paper consists of 
a necessary and sufficient proof of the existence of stable 
orbits, though only applicable, for the moment, to the 
lowest energy one for two different sets of closed trajec- 
tories. By stable we mean an orbit that, under the ef- 
fect of a (small) perturbation, tends naturally to go back 
to its initial configuration. Those perturbations would 
correspond, in this framework, to the instantaneous dif- 
ferences between emitted and absorbed radiation, once a 
stochastic background is included.. 

As a first step, in Sec. IIX1 a rigorous treatment is pro- 
vided only for the trivial case of circular orbits in the 
plane, for which we provide a necessary and sufficient 
condition. These orbits have, nevertheless, non-vanishing 
AM, and because of this we can presume they do not 
correspond to the ground state, but to the so-called p- 
orbitals. Another question is that those circular orbits 
are isotropic (around their axis of symmetry), while real 
p-orbitals are clearly not. This could be due to the fact 
that our simplistic model of the atomic nucleus assumes 
its isotropy, whereas the real hydrogen nucleus posses a 
clear break of isotropy inherent to the quark model. 

In Sec.[Xl an additional (infinite) set of stable solutions 
is provided to account for the whole (discrete) s-orbitals 
spectrum. These trajectories show vanishing (average) 
angular momentum. The main axis of each of these s- 
trajectories may randomly precess around the nucleus 
giving rise to the expected spherically symmetric quan- 
tum s-orbitals. 

Particularly in the case of our s-orbits, it is convenient 
to bear in mind that those trajectories are not stable, nor 
they are stationary, in the simplified problem of a scalar 
central potential. For them trajectories to be possible, 
we must add both the action of the Lorentz force on a 
point particle, as well as our higher order terms arising 
from the multipole orders we include in the calculation 
(our "stabilizing mechanism"). The identification of that 
last set of stable trajectories with vanishing AM, a set of 
"candidates" for the s-orbitals, is maybe one of the most 
important contributions in our present work. Indeed, to 
our knowledge, it is something not present in any previ- 
ous treatment on the subject (and of course not in the 
Bohr model!). 



G. Equivalence with QM orbitals 

Beyond the fact that this SED picture may give rise to 
stable orbits for the hydrogen atom, the equivalence with 
the well-known quantum orbitals has been established in 
the following ways: 



(i) Average angular momentum (AM) over the orbit: 
we have shown the existence of stable orbits whose aver- 
age AM either vanishes (for the s-orbits) or not (for our 
p-orbits). See Sees. H (TX] and EH 

(ii) Symmetry properties: as seen in Sec. IXVI C[ once 
phase averaging is done, rotational symmetry is recov- 
ered, three dimensional for the s-case and around an axis 
for the p-cases.. 

(iii) Parametrization of the discrete spectrum by a 
counterpart of the principal quantum number, obtained 
from an additional frequency matching condition. See 
Sees. IXIFl IXIH1I and IXIG1I 

H. Energy spectrum 

A remarkable approximation to the typical 1/n 2 , at 
least when we include only the principal oscillation, is 
also obtained. See Sees. IXIII and the Appendix, Sec. IB 51 
It is, anyway, a very preliminary result, given that the 
quantitative analysis is reduced to two dimensions, an 
isotropic model for the structure of the proton (very far 
from a realistic one) and the lack of a rigorous approach 
to the inclusion of the stochastic background. Some 
of these convenient extensions are discussed briefly in 
Sec. EX] 

I. A possible relation with Pauli statistics? 

This is an incidental question, but not devoid of in- 
terest at all. Basically, in our framework, stable orbits 
correspond to modes of oscillation. That correspondence 
is one to one, which suggest an interpretation in terms of 
a sort of "classical" Pauli principle. 

There is a very interesting interpretation of our results 
here. Stability of an orbit ( or the orbital generated from 
it) requires vibration of the nucleus in a certain mode. 
The presence of one of these modes can exclude (or not) 
the existence of more instances of the same mode, giving 
rise to a sort of Pauli principle that nevertheless we will 
only suggest here. For example, for the s-case everything 
seems to suggest that the only possible configuration for 
two electrons is the one in which their spins are anti- 
aligned, and that no third one can be included in the 
stable picture, for a given value of our principal number 
n s . The question would require, as said before, further 
attention. 



J. Extension to 3D 

A great deal of our conclusions are immediately gen- 
eralizable to 3D. For instance, we have justified three di- 
mensional rotational symmetry from our proposed stable 
s-orbits. However, in general, extension to the 3D prob- 
lem requires some extra calculations that for the moment 
we have not faced here. In particular, the effects of the 
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mechanism we expose here are clearly modulated by the 
classical counterpart of the LS interaction in QM, be- 
tween the electron and the nucleus, when the system is 
allowed to evolve in 3D. 



K. Additional conclusions 

With the aim of shedding some more light on our work, 
in Sec. lXViTl we have discussed the approach we adopted 
to include, in the dynamic equations of the system, the 
effect of rad/abs as a result of the interaction with the 
stochastic background. Certainly, our approach is not 
the most rigorous one, although it is enough so as to 
show how the equations can lead naturally to a set of 
privileged trajectories or attractors. 

Of course, all this only constitutes a preliminary ef- 
fort, and more detailed calculations about the specific 
predicted parameters of those stable orbits (or orbitals, 
in a wider sense) are awaiting to be done. Moreover, even 
for the case of elliptical orbits, a strict proof of stability 
is still lacking. The necessary elements are there, though, 
and at least in this particularly point, further work is on- 
going and will be presented elsewhere. A point that is 
convenient to stress is that our mechanism belongs to a 
"simplest one" class. Therefore, calculations based exclu- 
sively on our picture here do not guarantee for reasonable 
estimates, for instance of the sorp energy spectra. Noth- 
ing forbids, nevertheless, that further refinements of the 
model could adjust more to the experimental results. We 
have seen, already, that special relativity arises the need 
to assume a "discrete" charge (and magnetic moment) 
distribution for the proton, and space for much more so- 
phistication is clearly still open in this issue. 

As a last remark, David Regalado, student, has per- 
formed some calculations on the mean speed of the elec- 
tron in our proposed s-orbitals that agree quite well with 
the accepted values for the quantum case. This calcula- 
tions will be presented elsewhere. 



XIX. ON FUTURE DEVELOPMENTS 

The main problem that will occupy our time in the 
next future will consist of trying to extract, in a more 
rigorous way than what we have done, the well known 
dependent-spectrum from our framework. In this regard, 
the calculations presented here can only be considered as 
very preliminary estimates. 

For the moment, what we have proven is that the ex- 
istence of an infinite, discrete, spectrum is the natural 
consequence of our model, and we have even done some 
estimates that bear some resemblance with the real quan- 
tities. In any case, from here several main problems are 
to be addressed, each of them will presumably further 
correct that first estimate: 

(I) The extension of the model to 3D, that we have 
already treated tangentially in some sections. 



(II) The inclusion of anisotropies in the model: mainly, 
the evident anisotropy in the structure of the proton, and 
its lack of inversion symmetry. This correction should 
make our classical orbits (or sets of them, in the three 
dimensions, and their associated probability distribution 
once background noise is included) look more like the 
actual quantum orbitals. 

(III) On the other hand, we have not so far predicted 
any specifical probability distributions, which would in 
the last term define an "orbital". Indeed, a very inter- 
esting extension of that work would be to predict the 
density of probability that we would find considering a 
"broad band" (high frequency-cutoff) model for the fluc- 
tuations. Up to now, we have used a "low band" model 
where the value of the ZPF was averaged over an interval 
(we did this through the inclusion of the x _ hinctions, see 
Sec. lVIBj) . small enough not to mask the dynamics of the 
system as an electromechanical mechanism, so we could 
calculate the "privileged" trajectories or attractors. 

We are confident, therefore, we have offered at least 
the foundations for a more quantitative prediction in the 
near future, both for the spectrum and the actual shape 
of the orbitals. 



XX. LAST PERSONAL REMARKS 

As a concluding remark, we must say that this whole 
work was born from an already old idea by the author. 
Curiously enough, in spite of my already old interest in 
SED, I was only partially aware of some of Cole's work. 
How close his approach to the issue of the stability of 
the hydrogen atom did not appear clear to us until much 
of the work was already in progress. It was especially 
delightful to realize how well our main idea (that mech- 
anism of energy exchange with an "internal" degree of 
freedom) fitted within the framework that Cole and Zou 
had settled in their already impressive amount of work. 
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Appendix A: Additional calc. on the IMSO 

We include here several additional calculations, regard- 
ing our idealized situation, with no necessary mention to 
the actual hydrogen atom. Therefore, notation here is 
in coherence with that of the first part of the paper: we 
always use R or b and W( ) because they are more general 
and perfectly defined variables (in contrast to Rq, Wo> . . ., 
which have different meaning for s or p-orbits) . 

Though not relevant enough to be included in the main 
body of the paper, they may be useful, however, to in- 
terpret further results. 



1. E-correction for p-orbits 

The aim of this section was to do a preliminary explo- 
ration: it seems, up to now, that the high order terms 
under our focus do not produce any correction... on av- 
erage over a whole circular (closed) orbit. This seems 
reasonable as all the new terms stem from the Lorentz 
force, for which "work" always vanishes (it is normal to 
velocity). This question should be treated more carefully 
in future versions of the paper. 

We want to evaluate what shifts in energy are intro- 
duced by our IMSO coupling. It is important to state 
again here that this new correction is of a completely 
different nature from the well known quantum LS. We 
have a coupling between (oscillating) magnitudes of clas- 
sical angular momenta, not just the one between their 
relative orientations in space. 

This said, we will confront the fact that (in some ways 
"as expected"), our correction vanishes over a whole 
cyclic orbit. This is a natural consequence of our cal- 
culations. Nevertheless, it would be interesting to think 
how this could be reformulated in terms of work, taking 
into account the special nature of the Lorentz forces. 

Specifically, we insist on the fact that a moving charge 
in an electrostatic field will "see" a magnetic field, due to 
Lorentz covariance. This magnetic field produces a force 
that couples to the charge. In a purely classical picture, 
we have seen we need to go up to a second order moment 
to be able to couple this field to the AM with which that 
charge is rotating around its axis of symmetry. 

We have, therefore, a classical coupling between the 
magnitude of the orbital angular momentum and of what 
we have called, with the aim of accommodating our rea- 
soning to a physical situation that could resemble a re- 
alistic hydrogen atom, "residual spin". Hence, after all 
what we have is a new correction term to the classical 
Hamiltonian (in the same way that one adds new terms 



to the quantum one). We focus now only on this new 
energy correction. 

On the other hand, and based on the formal equiva- 
lence of the equations exposed in Sec. IV Bl we can con- 
sider that the following calculations are valid for either 
of the treated situations, the first one comprising two 
systems of charges, a heavy pointlike one and a lighter 
extended one, the second being a first model of the hy- 
drogen atom. 

As a first approach to characterize quantitatively this 
new energy correction, we begin by equating, classically, 
the (radial) force with the centrifugal term. We will first 
focus on the order of the correction (then more detailed 
conclusions can be obtained by using the electron, proton 
mass, charge, magnetic moment, etc). 

For a perfectly circular orbit, we have, from (|72p with 

Vr = and the perfect compensation (as an expectation 
value) of radiative loss and background pressure, 



F P + fW + 



The terms F e and Fm •* are already present at any clas- 
sical treatment of the hydrogen atom, and represent, re- 
spectively, the electrostatic and the Lorentz (reciprocal, 
comment on this!) forces. The terms that give rise to 
F^ are the interesting ones. 



(°)\2 



Ro 



(Al) 



a. Kinetic E-correction 



The kinetic energy T = ^ m (v^) 2 , and we can com- 



pute each correction by: 



ACT), « F t , 



(A2) 



Specifically: 



^ oc Rj b AT e oc R- r \, (A3) 
AT^^R- r \, (A4) 



F£> oc qB 1 vi o) oc Rj b 
F m oc — — oc R~J b => AT(f ) cc R- r l (A5) 



R 



orb 



where we use the sign A all the time because a complete 
evaluation of the kinetic energy would have to include 
the term in w 2 . From the perturbative point of view, we 
are only interested in the former ones. 



b. Potential E-correction 



We need to compute AR or b as a result of F^, and 
then do, to first order, 



AV oc 



1 



1 



Rorb Rorb + AR or b 



(A6) 
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This accounts for solving a polynomial equation. We 
multiply by R orb + AR orb on the right, 



{Roi 



AR 



AV oc 



AR. 



orb 



{Rorb 

and therefore 



Rorb 

Rorb ■ AR orb ) AV oc AR orb , 



AV oc AR 



but we also had, from (|A5|) . AT^ oc R a ^ b - Now, in first 
order to perturbation (we consider the other terms of the 

kinetic energy T e and T$ remain fixed), again selecting 
the leading terms, 



■b- 



(A7) 
(A8) 

(A9) 



ATI*) 



J orb 



[ {Rorb + AR or b) 2 ~ R or b ] 



oc -m u? orh [ 2R orb AR orb 

OC m Rorb ARorb, 



and therefore, 



AR orb oc R„ 



(A10) 



(All) 



so we can neglect the correction to the potential energy 
AV oc R- r 6 b . 



c. Order of the full correction 

Our first idea was to compare our estimate with the 
value of {ipGs\ L ■ S \i^gs), the first correction in the 
perturbation expansion of the energy predicted by non- 
relativistic QM for the ground state \tpGs) of the system. 
Now, from [32J (page 88), we have, for a central potential, 
the leading term with L ■ S has the magnitude 



1 dV 
— oc R 

r ar 



3 



Ron 



VkR: 



(A12) 



for a Coulombian potential. It does agree, therefore, with 
the order of our leading correction ATm (!) in (|A5I) . 
Nevertheless, we must stress our correction is of a com- 
pletely different nature, coupling oscillations on the mag- 
nitude of a classical AM vector, not on its orientation as 
the quantum LS does (see Sec. Ill D|) . 

Then, what is the interest of this comparison? In nu- 
clear physics, phenomenological (ad hoc) terms with the 
LS form are introduced in order to explain the negative 
binding energies (therefore, the stability) of most nuclei. 
This is a very interesting question that we should treat 
elsewhere. 

There was, anyway, a triviality waiting to be acknowl- 
edged here: for a perfectly circular orbit, we simply have 
AT„r = 0. But still we were interested in that calcula- 
tion, due to, for an elliptical orbit, it still remains valid in 



the points where the radial velocity vanishes. In the rest 
of the points, our correction to the kinetic energy needs 
to include the velocity component in the radial direction, 



Vi- 



ta) 



From (|72p we have 



\4 o) \ocR 



-4 
b 



\v^\ dt<xR- or \ 



ATocb(°)| 2 oc R-* 



(A13) 



so therefore we can neglect that contribution and simply 
extend our former analysis to any possible orbit (we have 
only strictly proved the existence of circular ones) . 

All this analysis justifies us to stick just to the leading 
term, and give a complete estimate, using (j37]): 



AT (R) 



Ro 



2 

3/tQ2 

Z ^orb 
3/tQ 2 

Z "-orb 



p(R) 



[ 



J {o) 



(A14) 



where in the last line we have used ©. To compute the 
correction, we would average the quantity ATm on all 
the points of the orbit. We do not tackle this calcula- 
tion because we have not yet found those stable elliptical 
orbits. 



Appendix B: Additional calculations on the H- model 

Here we provide additional calculations on the model, 
once this is completely particularized to account for the 
hydrogen properties. Whenever the results here are used 
in the body of the paper, it is always properly referenced. 
Notation here corresponds to that of the last part of the 
paper: for convenience, and because there is no possible 
ambiguity, we use Ro,loo,Ri, instead of R^ , Rq , . . .. 



1. Energies without 2nd-Osc 

a. Energy of a stationary p-orbit 

A circular stationary trajectory A p is given by the pair 

of parameters ujq = uj( ) arL d Rq = R or b- Assuming vi ^ — 
0, we can accept the following expression, from (|B3[) . 



R'n 



1 



4ire m e u) 2 ' 



and therefore 



R oc (wo) 3 , 



(Bl) 



(B2) 



which means that, in general, if ujq decreases, Rq will in- 
crease, which is intuitive for a higher energy orbit. After 
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the inclusion of the secondary oscillation the situation 
will be, however, slightly more subtle than that, with i? 
keeping an almost constant (and perhaps even decreas- 
ing) value, that nevertheless will not be associated with 
a decrease in energy. Now, for a circular stationary orbit, 

v[°^ = u)t \Ro, and we have 



(w (o) ) 2 
rn P _ — — — = m e u> Ro 



Ro - e ~ u " u 47reoJ^' 
and, therefore, in terms of energy, 



(B3) 



1 ( (o)\2 

-m e (vl ') 



1 

rn e l 

1 q 2 



2 AneoRo 



~qV(Ro) = -~E vot , 



(B4) 



with V(R or b) the potential created by a charge with value 
q > at r = (ri = — R or i,xi), and therefore, with, 



E = E, 



kin 



E, 



pot 



pot 



Ep t — —E, 



pot J 



(B5) 



where both E k i n and E pot stay constant along the orbit. 
On the other hand, we can also reverse (|B4[) . obtaining 



E, 



pot 



and finally, 



E 
AE 



-AE ki 



-E, 



pot ■ 



(B6) 

(B7) 
(B8) 



b. Energy of a stationary s-orbit 

On the other hand, an s-type trajectory A s is for in- 
stance determined by the pair of parameters v max = 

maxA s Vr and R max — maxA, Rorb = Ro (this last 
equivalence only under the approximation of harmonic 
dynamics both in the radial and tangential directions). 
Now, we can equate kinetic and potential energies (for 
this second one we adopt a convenient absolute reference 
so it vanishes when it is minimum) at the extremal points 
(when the velocity is maximum and when it is zero): 

E Hn {R = 0) + E pot (R = 0) = E pot (R = Ro), (B9) 

and (with I z the z-moment of inertia of the nucleus), 

E kin {R = 0) = im e (w i? ) 2 + ^I z cj 2 , (BIO) 



and 



E pot (R = R ) = -qV{R = R )-E pot (R = 0) 



47re i?o 



— E pot (R — 0), 



(Bll) 



and finally, 

E po t(R — Ro) = Ekin{R = 0) + E pot (R = 0), 

I 

Ekin(R = 0) = E pot (R — Ro) — E pot (R = 0), 



(B12) 



(B13) 

but now, disregarding the contribution to the kinetic en- 
ergy coming from w, 



1 q 2 
-m e (w i?o) 2 = -: 5-, 

2 47reo.Ro 



(B14) 



and indeed, if no other contribution to the energy is con- 
sidered (no oscillation in v\ or uj are present), the aver- 
age energy is 



Ekin(R — 0) — —E P ot(R = Ro) 
E = -E kin (R = 0), 
1 
2 



AE — —AEkin — 77 E pot ■ 



(B15) 
(B16) 

(B17) 



Of course, the former calculation must assume that the 
potential energy is bound in the origin r = 0, somehow. 
This is however compatible with one of the main ideas 
behind this paper: the introduction of the inner structure 
of the nucleus in the atomic model. If we adopt now the 
harmonic oscillator approximation, with R = Rq cos(cjo^) 
(where from ujq does not stand for a rotation in respect 
to RFq but can still be defined as an angular frequency) , 
then v r \ max — RoLlIo- Also working with (|B14j) . 



^m e (uj Q R Q ) 2 

Rq 



47re i?o ' 
2ireom e Lo 2 ' 



(B18) 
(B19) 



and again to arrive at a similar condition as in the circular 
case, this time for the parameter Rq, 



R 2 (x(lu ) 3. 



Potential E-gap between adjacent orbits 



(B20) 



We recall once more that, either for the p or s case, 
and with no rad/abs, our model admits an infinite set 
of stable trajectories, therefore an infinite set of pairs 
U( ),R rb and v max ,R max . That continuous set will be 
reduced to a discrete one by means of the imposition 
of that balance. Let us now consider (feasible) excited 
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states. The main idea here is that, either for the s or p 

cases, 



AE oc AE 



pOt-! 



(B21) 



and moreover, if R or b {R ma x) increases, AE > because 
AE pot > as E pot < 0. 

See where we are going? The same spectrum for s and 
p-type trajectories demands R or b ^ Rmax for a given 
value of the principal quantum number n, but this is 
no problem. Let us, for the moment, just work with a 
generic parameter R n (corresponding to an excited state 
with quantum number n). It is interesting to quantify 
this increment between the trajectories n and n + 1 as 



AE oc AE 



pot 



AlTEr 



1 



Rn 



1 

R„ 



(B22) 



and 



r m 2 ui( R ) l 

L * K n J J max 

M 2 u>r,o 

(o)| 



3/i e Q 2 _ (o) 

Kg 

3/t e Q 2 



7? 4 



(B25) 
(B26) 



naturally with v r \ m ax — Ri^i- The two former equa- 
tions allow us to eliminate oj r $ and produce a relation 
between Ri and u±, which is what we were looking for: 



3fi e Q 2 _ 9//gQ 2 Q 2 
* ~ M 2 R* 



(*) 



m e Rfj 



3ft e ^J Q 2 Q i 2 * ) /m e M : 



(B27) 
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2. Determining the parameters of 2nd-Osc 



2nd-Osc for s-orbits 



Aside from the harmonic oscillator approximation, we 
neglect the change in the electrostatic and Lorentz parts 
of the potential as a result of the change in the ra- 
dial coordinate R = \R or b\- Those electrostatic and 
Lorentz terms contribute to the mean centripetal acceler- 
ation, but nothing else, at least under our approximation. 
Therefore, only the term corresponding to our high order 
model enters in the following calculations. 



a. 2nd-Osc for p-orbits 

For circular orbits, the equations of interest are ([72]) 
and (fT4")l . Considering that, as already stated, the sec- 
ondary oscillation (2nd-Osc) can be expressed assuming 
the radial direction as AR = R — Rq = R\ coswii, this 
oscillation univocally determines an associated one for 
ur = u)r R \ (see footnote [HI): 



lo r = -uj + w R<0 cos(wit), 



(B23) 



where uj r .q is determined from ljq, RoUJi, Ri. Now, we 
can apply the following two conditions, evaluating the 
point of maximum acceleration, measured as a variation 
of the radial velocity (we are already eliminating the cen- 
tripetal contribution, associated to ujq in the former ex- 
pression): 



m e vi o) 






max 



3Ae 



m e RiUJ 1 



3fi e Q 2 



(B24) 



For circular orbits, the equations of interest are ([72]) 
and ([73"1) . Considering that, as already stated, the sec- 
ondary oscillation can be expressed for the tangential di- 
rection as R or b -yi = R\ cos uj\t. As it already happened 
for circular orbit, this oscillation univocally determines 
another one, in this case associated with t>r°': 

= uj R o cos(w t) + Av(°\ 
Avl o) = Atv,ocos(u;ii), 



(B29) 
(B30) 



where Aiv,o is determined from toi , R\ . 

Now, we can apply the following two conditions, eval- 
uating the point of maximum acceleration, measured as 
a variation of the radial velocity (again we are already 
eliminating the main, centripetal contribution): 



•0) 

m e v; 



max 
1 



and 



A(m^ o) ) 

m e ujiAv r fi 
m e Av r o 



qB x Avi o) 
qBi Av rt0 , 



qBt v { t o) 

qB 1 Riui, 
qBi Rl 



(B31) 



(B32) 
(B33) 



Now, the two equations allow us to eliminate Au r ,o and 
produce a relation between R\ and u)\ , which is what we 
were looking for, 



ui = 



ggi 

m 

qve 

R 3 ' 



q b x 



2 ~ 2 

q k 



(B34) 
(B35) 
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3. Approach to the orbital radius spectrum 

a. p-orbits 

We know 



and so 



3 1 2 1 111 

R oc — , R oc — , R oc — , —ocn$. 
n„ „3 3 iio 



p 



'p 



(B45) 



m^Ro 



Ro 



and therefore 



R n 



(B46) 



UJq = LJl, 



(B36) 



and using (|B28[) as well as our external (it is a supple- 
mentary assumption) "matching condition" ojq = — wi, 



1 3Aa\/gFA 



4ire rt 



,(*) 



5 4-n-ep 



so therefore 



(B37) 

(B38) 
(B39) 



(B40) 



c. Comments 

It is surprising that Rq is actually a decreasing func- 
tion of n p ,n s . Nevertheless, what we do here is simply 
provide a feasible stable orbit, from which a realistic prob- 
ability distribution ( extended, in principle, to the whole 
space) could be generated once the stochastic background 
is included in the picture. Actually, the inner regions or 
the quantum s and p-orbitals do tend to decrease! Be- 
sides, in a real atom repulsion between electrons can well 
explain that the actual mean distances in the orbitals 
increase with the principal quantum number. 



4. Modified power balance 

So far, we have seen that given Rq, luq and u>i are de- 
termined (the second through our "matching" condition) . 
Therefore, the only free parameter is R\. Imposing the 
radiated/absrobed power balance will fix it. 



and so 



W CX — 3 OC Up . 
Tin 



b. s-orbits 



(B41) 



Under the harmonic oscillator approximation, the 
spectrum of the pendulum orbit maximum radius fol- 
lows the same rule as for circular orbits. Indeed, from 
the stationarity equation, where we disregard higher or- 
der terms: 



m e w%Ra 



AneoRo 



qti e 

Rq" 



1 

2n, 



-wi, 



(B42) 



and using (1B28[) and the matching condition u>q = ^-toi, 



1 gAe 
2n p Rq 3 



Ro = 



R'a 



4vreoi?§ ' 



(B43) 
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a. p-orbits: first approach 



Firs, we notice that ([78)) and (|80[) provided instant av- 
erages, rather than averages over the whole orbit, which 
is important to bear in mind for the following calculation 
(otherwise, due to ujq < u>±, some extra factors would be 
necessary). Now, recovering (|78|) and (|80|) . we have, 



(Prad) 



„2, ,4 d2 

67reoc 3 

\2tT€qC : 



127re c 3 ' 



[2w*R 2 + R^tRj] 



as well as 
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67reoTO e c 3 
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q 2 hujl 
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(B47) 



(B48) 
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Now, applying P rad = P ab and u>i = n p u> , 



s-orbits 



2utR 2 



III' P. 



2u*R 2 + n%u\R\ = — ■ [2ul + n 3 p uj 



3 i] , 



2w i?o + nfooRl = [2 + np] 



and 



R 2 



— ■(2 + nl)- 2lo R 2 



(B49) 
(B50) 
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(B52) 



and now we can do a preliminary calculation regarding 
only the highest order terms, and in the limit where n p — ► 
oo, taking into account that Rq oc —jjs, 



R\ 



1 



R 2 -> 



1 ft 



— < + 2uj q R 2 



(B53) 
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n p w Tie 
and therefore i?i — > when n p — > 00. 

6. p-orbits: detailed calculation 

A rigorous calculation of absorbed power P a ;, for the p- 
case implies taking into account the following expression 
for the trajectory: 

x = Ro [ 1 + Ri cos(wit) ] cos(wot), (B55) 
y = i?o [ 1 + Ri cos(wit) ] cos(wo^)7 (B56) 

and now we would have to write the radiation damp- 
ing equation for each of this components. An easy path 
can be developing the cosine products in an addition of 
cosines: we would have the sum of three harmonic oscil- 
lators for each x or y dimension. 

(complete) 

On the other hand, the former expression for radi- 
ated power can be justified starting from the following 
general formula for instantaneous loss due to radiation 
( Abraham-Lorentz , non-relativistic) : 



P' 



Hoq 2 a 2 

QlTC 



where a is the acceleration modulus, 
orbits we would have 



K (o) ) 2 



•,(0) 



R, 



orb 



(B57) 

Now, for our p- 
(B58) 



from where we can see that, when doing the square and 
integrating along a whole closed orbit, cross terms clearly 
vanish, and we have just the addition of two independent 
contributions. 

(complete) 



For s-orbits, no departure from Puthoff's initial calcu- 
lation, except for the inclusion of a factor, is necessary, 
with 



x = Rocos(uiot), 
y = i?i cos(wit), 



(B59) 
(B60) 



and this time we have the typical case of two harmonic os- 
cillators in quadrature. The radiation damping equation 
for on oscillator in one dimension is well known, and the 
calculation can follow Puthoff's steps making the neces- 
sary distinctions. We have, at the end, two contributions, 
one coming from the radial velocity and the other coming 
from the tangential one. The calculation is similar to the 
p-case, but some factors change, as now the frequency 
matching condition is uj\ — 2n s u)Q, to keep zero average 
AM. Applying P rad = P ab and oj x = n p iv , 



^Rl + ufRl 



oj^Rl + IGnj^Rl 
uqRo + 16n^ujQR 2 
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h 
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lo a + 8n>g , (B62) 
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and 



1 



( 1 + 8nl ) - 2lu q R 2 
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5. 2nd-Osc in the E-spectrum? 



In IXII C[ we have commented on why we would not 
include the 2nd-Osc into the energy spectrum. Here we 
nevertheless do some extra calculations to have better 
elements of judgement. Before proceeding, we recall that 
once the secondary oscillation is included in the picture, 
there are (at least) two main contributions to consider: 
(i) a kinetic energy coming from the movement of the 
electron as a point particle (contributions to the potential 
energy vanish on average for this movement), and (ii) a 
kinetic contribution arising from oscillations in uj, i.e., 
the nucleus (proton) rotating as a whole. As usual, the 
sum of the two terms is constant, and at the extreme 
points of this oscillation, always one of them vanishes, 
simplifying the calculation. 

In particular, following our calculations in Appendix 
IB 41 eq. (|B28I) , we could work on the power balance equa- 
tions of the p-orbits, estimating the (kinetic) energy of 
that secondary oscillation. From (|B49|) , 



2uqRI 
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(B65) 
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and now we can isolate 

Z Z LU-y 

(B66) 

E 2 (n p ) = ]-m e ojfRj = -— [h(2 + n^)u Q - m e ulR%\ , 
Z Z n p 

(B67) 

where the only interesting term is the one that grows like 
npijQ. This is not a cause of concern, however, because 
luq << m e ujQRQ in general, and we should also bear in 
mind that all this is no more than a very gross approxi- 
mation. Besides, at higher frequencies higher order terms 
in the multipolar development must take action, surely 
with the consequence that a high-frequency cut-off can be 
effectively introduced in the model. That cut-off would 
eliminate the problem that the term in uj\ poses for our 
prediction of the spectrum. 

Also, with m e ujQRl « H (now Puthoff's condition is 
not exactly fulfilled for the orbital movement, but we 
may still be able to approximate), 

E 2 (n p ) « \^ m ^\R\ « ^m e uj^Rl, (B68) 
where the variables on the right hand side, Rq, u>q, etc 



would also depend on n p . Also, from (|B61|) . we do the 
same for the s-orbits, arriving to a similar result, 

E 2 (n s ) m )-m & u\R\ pb ^m e Lj%R 2 , (B69) 

but, as seen in Appendix [B] either for circular or pendu- 
lum orbits, this means precisely 

E 2 ~ E kin .i = —Ei, (B70) 

and therefore here, surprisingly only in principle (read 
below), with n = n p ,n s , the gap between adjacent 
states obtained from the primary oscillation would sim- 
ply get cancelled out: AE 2 (n) = E 2 (n + 1 ) — E 2 (n) = 
-AE 1 (n,n + 1). 

The interpretation of this last result is still 
an open question for us. Our guess is that 
the secondary oscillation must be regarded as 
as a phenomenon that is "transparent" to QM: 
there needs to be, certainly, some secondary wave-packet 
emission/absorption associated to an atomic transition, 
though at a higher range of frequency ( such high frequency 
perhaps explains why they are not detected, remaining 
completely transparent to ordinary QM). 
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